one  FILE  COPY  ADA  1 12  7  ol: 


AFAPL-TR-78-81 
Volume  II 


FOREIGN  OBJECT  IMPACT  DESIGN  CRITERIA 


Albert  F.  Storace 


General  Electric  Company 
Aircraft  Engine  Group 
Cincinnati,  Ohio  45215 


February  1982 

Interim  Report  for  Period  31  December  1977  -  31  December  1978 


Approved  for  public  release;  distribution  unlimited. 


AERO  PROPULSION  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT -PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


DT1C 


elec 

MAR  2  2 


TE 

1982 


B 


82 


03  22 


001 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  In  connection  with  a  definitely  related  Government  procure¬ 
ment  operation,  the  United  States  Government  thereby  Incurs  no  responsibility 
nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  In  any  way  supplied  the  said  drawings,  specifications 
or  other  data,  Is  not  to  be  regarded  by  Implication  or  otherwise  as  In  any 
manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  Invention 
that  may  In  any  way  be  related  thereto. 

This  technical  report  has  been  reviewed  and  Is  approved  for  publication. 


Project  Engineer  Technical  Area  Manager, 

Propulsion  Mechanical  Design 


FOR  THE  COMMANDER: 


.  SHIPMAN,  MAJ,  USAF 
Chief,  Components  Branch 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  If  the  addressee  is  no  longer  employed  by  your  organization,  please  notify 
AFWAi/PCrrc  WPAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list. 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obi Igatlons,  or  notice  on  a  specific 
document. 

AIP  roncc/s«7a0/12  M,rch  1912  —  70 


_ UNCLASSIFIED _ 

SSCURITV  CLASSIFICATION  OF  THIS  RAPS  (Wltmm  Ornm  toMnQ _ 

REPORT  DOCUMENTATION  PAGE 

I.  RSRORT  NUMOKR  II.  OOVT  ACCSSMON  Moi 

AFAPL-TR- 78-81,  VOL  II  JlP  ~A  1 1  '2-  1( 

4.  TITUS  (tmT  SumnUrn) 

FOREIGN  OBJECT  IMPACT  DESIGN  CRITERIA 

7.  AUTHOR?*) 

Albert  F.  S Corace 


uao  MtntucciOMs 

BBXHUt  COMPLmWO  FOWI 
:ciPiiNri  cataloo  nunscr 


*■  tym  or  rcrort  •  rcieoo  covsrco 

Interim  Technical  Report 
31  Dec.  1977  -  31  Dec.  1978 

«.  rieroMNNS  om.  rcrort  nuhrsr 
R78AEG180 


f.  Contract  or 


'  noBRUWfJ 


F33615-77-C-5221 


renroMMiNO  organization  nans  amo  aoorsss 
General  Electric  Company 
Aircraft  Engine  Group 

Cincinnati,  Ohio  45215 _ 

II.  CONTROLLING  OFFICE  NANS  ANO  ADORS1S 

Aero  Propulsion  Laboratory  (AFWAI/POTC) 

Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  Air  Force  Base.  Ohio  45433 
TT  monitoring  aocncy  namc  S  aoorssv"  SSESt STeSniSS  o««»; 


rrooram  element.  froject.  tamT 
AIWA  *  ww  unit  numosrs 


3066-12-33 


It.  REPORT  OATS 


February  1982 

•a.  MuSeiii  or  r Sen 

72 _ 

is.  ascuRiTY  class.  (mi  m*  mmhj 


IE  distribution  statement  (mi  wm  Rapao 


Unclassified 


fi  cation/  oown  grading 


Approved  for  public  release;  distribution  unlimited. 


I  IT.  DISTRIBUTION  ar ATUfleNT  (ml  »m  Allwl  nM  W  BIm*  19,  ll  mUmmmmt  *mm  Kmmmt) 


is.  supplementary  worca 


The  computer  program  contained  in  this  technical  report  is 
theoretical  and  in  no  way  reflects  any  Air  Force  owned 
software  programs. 


It.  KCY  WORDS  /CnmRnm  mm  rmrmttm  im>  If  muMr  i 

Transient  Structural  Response 
Foreign  Object  Loading 
FOD  Design  Criteria 
Damage-Tolerant  Design 


Fan  /Compressor  Blad-*  » 
Material  Transient  -»  » 

Structural  Failure  ‘..r^...*  > 

Airfoil  Structural  Modeling 


»  Alltl/lfT  (Cmttmm  m  wyw—  NMt  If  mhi»oi>  n  My  ll—l  mmbm) 

~The  program  objective  is  to  establish  specific  design  criteria  and 
provide  the  analytical  design  tools  to  assess  and  improve  the  foreign 
object  damage  tolerance  of  turbine  engine  fan  and  compressor  blading.  -> 

v  This  program  will  aid  in  the  design  of  more  efficient  damage-tolerant 
blading  through  the  replacement  of  trial  and  error  FOD  test  and  evaluation 
practices  with  systematic  transient  structural  analysis  methods,  test  pro¬ 
cedures  and  design  criteria.  — _ _ _ _ 


SDIUON  or  I  NOV  M  IS  OeSOLtTI 


_ UNCLASSIFIED 

ICCURITT  CLASSIFICATION  Of  TNII  RAGS  i 


UNCLASSIFIED 

iccujjTY  d.*Miyg«rtow  or  thi*  rAoymum 


design  system  structure  was  collated  and  Initial  versions  of  the 
preliminary  end  final  design  trsnsient  structural  response  computer  progress 
were  developed.  In  addition,  the  interface  between  the  bird  loading  aodel 
and  the  final  design  response  aodel  was  coapleted.  Three  first  stage  blades 
were  selected  for  modeling  and  testing.  Tests  and  analyses  to  be  used  to 
generate  local  and  gross  structural  dsasge  design  data  ara  in  progress  and 
planning  was  coapleted  for  structural  element  impact  tests.  Composite  and 
aetallic  specimens  for  local  and  gross  damage  testing  were  designed  and  most 
of  these  specimens  have  been  fabricated.  A  parametric  matrix  to  guide  the 
analyses  and  validation  testing  was  completed  and  FDD  design  criteria  goals 
were  formulated  for  consideration  in  the  program. 


UNCLASSIFIED 


mcukitv  ck«Miric«rioM  or  tmi*  i 


PREFACE 


This  report  covers  work  performed  by  the  Aircraft  Engine  Group  at 
General  Electric,  Evendale,  Ohio.  This  report  covers  work  completed 
during  the  period  1  January  1978  through  31  December  1978.  The  contract, 
F33615-77-C-5221  was  sponsored  by  the  Aero  Propulsion  Laboratory,  Air 
Force  Systems  Command,  Wright-Patterson  AFB,  Ohio  45433,  under  Program 
Element  62203F,  Project  3066,  Task  Area  12,  Work  Unit  33  and  under  the 
direction  of  Sandra  K.  Drake  (AFWAL/POTC) ,  Project  Engineer.  Albert  F. 
Storace  was  technically  responsible  for  the  work  and  preparation  of 
this  report  in  fulfillment  of  the  aforementioned  contract. 


Accession  For 

KTTS  C^A&I 
MIC  TiP, 
Unamoutvei 
Ja  ji  . :  i i  ;  n 


□ 

□ 


By - 

1  Distribution/ 


Availability  Codas 
(Avail  and/or  — 


Dist 


Special 


J 


TABLE  OF  CONTENTS 


Section 

1.0 

2.0 

3.0 

4.0 


5.0 


6.0 


INTRODUCTION 


SUMMARY 


2 


TASK  I  -  DESIGN  SYSTEM  STRUCTURE  3 

TASK  II  -  TRANSIENT  RESPONSE  ANALYSIS  MODELS  6 

4.1  Introduction  6 

4.2  First-Level  Response  Model  (NOSAPM)  7 

4.3  Second -Level  Response  Model  (COMET-BLADE)  11 

4.4  Interface  With  the  Loading  Model  11 

TASK  III  -  IMPACT  LOADING  MODEL  15 

5 . 1  Introduction  15 

5.2  Subtask  A  -  Simulate  Bird  and  Ice  Impact  Parameter  15 

5.3  Subtask  B  -  Loading  Model  Sensitivity  Study  15 

5.4  Subtask  C  -  Formulation/Interfacing  of  Loading  Models  With  17 

Structural  Models 

5.4.1  Introduction  17 

5.4.2  The  Slicing  Model  18 

5. 4. 2.1  Slicing  Model  Development  for  Birds  18 

5. 4. 2. 2  Slice  Geometry  Parametric  Study  22 

5.4.3  Fluid  Flow  Model  of  Impact  22 

5. 4. 3.1  Surface  Singularity  Method  -  Superposition  23 
of  Onset  Flow  and  Surface  Source  Distribution 

5.4. 3.2  The  Body  Surface  Approximation  23 

5. 4. 3. 3  Solution  Procedure  23 

5. 4. 3. 4  Onset  Flow  Modeling  24 

5.4.4  Structure  of  the  Loading  Model  and  Interface  With  24 

Response  Models 

5.4.5  Summary  26 

5.5  Subtask  D  -  Experimental  Verification  Loading  Model  26 

TASK  IV  -  MATERIAL  RESPONSE  AND  FAILURE  CRITERIA  28 

6.1  Subtask  A  -  Gross  Structural  Damage  Properties  28 

6.1.1  Introduction  28 

6.1.2  Real  Blade  Strain  Rate  29 

6.1.3  Material  Property  Data  31 

6.2  Subtask  B  -  Local  Leading  Edge  Damage  41 

6.2.1  Introduction  41 

6.2.2  Overall  Approach  43 

6.2.3  Metal  Specimen  Tests  43 


v 


TABLE  OF  CONTENTS  (Concluded) 


Section  Page 

6.2.4  Composite  Specimen  Tests  44 

6. 2. 4.1  Laboratory  Size  Composite  Speclaan  Tests  44 

6. 2. 4. 2  Full-Scale  Composite  Blade  Tests  44 

6.3  Subtask  C  -  Material  Screening  Tests  44 

6.3.1  Introduction  44 

6.3.2  Materials  to  be  Investigated  45 

6.3.3  Screening  Tests  45 

6. 3.4.1  Ballistic  Limit  Testing  45 

6. 3. 3. 2  Maximum  Deformation  Testing  46 

6. 3. 3. 3  Gross  Structural  Damage  Testing  46 

6. 3. 3. 4  Fatigue  Testing  46 

7.0  TASK  V  -  PARAMETRIC  ANALYSIS  47 

7.1  Introduction  47 

7.2  Blades  Selected  for  Analysis  47 

7.3  Parametric  Matrix  47 

7.4  Modeling  of  the  Structural  Specimens  51 

8.0  TASK  VI  -  STRUCTURAL  ELEMENT  TESTS  56 

8 . 1  Introduction  56 

8.2  Composite  Specimen  Tests  56 

8.2.1  Laboratory  Size  Composite  Specimen  Tests  56 

8.2.2  Full-Scale  Composite  Blade  Tests  58 

8.3  Metal  Specimen  Tests  58 

9.0  TASK  VII  -  ERROR  BAND  ANALYSIS  61 

9.1  Introduction  61 

9.2  Properties  Requiring  an  Error  Band  Estimation  61 

10.0  TASK  VIII  -  FOREIGN  OBJECT  IMPACT  DESIGN  CRITERIA  63 

10.1  Introduction  63 

10.2  Foreign  Object  Impact  Design  Criteria  Goals  63 

11.0  REFERENCES  64 


vi 


LIST  OF  ILLUSTRATIONS 


Figure  No.  Page 

1.  Program/Tasks  Flow  -  Foreign  Object  Impact  Design  Data.  4 

2.  Foreign  Object  Impact  Design  Criteria  Program  Functional  5 

Organization. 

3(a).  Level  1  Response  Model  Formulated  With  3-D  Finite  Elements.  8 

3(b).  3-D  Finite  Element  Used  In  NOSAPM.  9 

4.  Plate  Modeled  From  NOSAPM  With  20-Noded  Elements  Rotating  at  10 

5000  RPM  and  Loaded  With  a  Transient  Load. 

5.  Level  2  Transient  Response  Model  Formulated  From  a  System  of  12 

Plate  Elements. 

6.  Link-up  Between  the  Loading  Model  Subroutine  and  the  Executive  14 

or  Main  Response  Model *Program. 

7.  Loading  Model  Sensitivity  Study  9-Element  20-Noded  Brick  Model  16 

of  Cantilever  Plate  Using  Symmetry. 

8.  Bird-Blade/Interaction  Geometry.  19 

9.  Impact  Test  Setup  for  Real  Blades.  30 

10.  J79  Stage  1  Compressor  Blade  Strain  Gage  Locations  for  Maximum  31 

Strain  Rate  Determination. 

11.  Strain  and  Strain  Rate  Time  Histories  for  Gage  No.  1  (at  Root  32 

of  Airfoil  on  Pressure  Side) . 

12.  Strain  and  Strain  Rate  Time  Histories  for  Gage  No.  3.  33 

13.  Strain  in  the  Time  and  Frequency  Domain  Resulting  from  34 

Hammer  Blow  for  Gage  No.  4  J79  Blade,  Test  Date  10/16/78. 

14.  Tensile  Specimen  for  Metals.  37 

15.  Split  Hopkinson  Bar  Specimen.  38 

16.  Modified  ITTRI  Specimen.  39 

17.  Modified  10°  Off-Axis  Shear  Specimen.  40 

18.  Through-the-Thlckness  Test  Specimen.  42 

19.  Comparison  of  20-Noded  Finite  Element  Solution  With  Beam  53 

Theory. 

20.  Comparison  of  16-Noded  Finite  Element  Solution  With  20-Noded  54 

Element  Solution. 

21.  Task  VI  B/Al  Composite  Test  Specimen  Configurations.  57 

vii 


LIST  OF  TABLES 


Table  No .  Page 

1.  Blade  Calculation  Input  Data  Sheet  £or  3 -Ounce  Bird.  20 

2.  Blade  Calculation  Input  Data  Sheet  for  1.5-Pound  Bird.  21 

3.  Characterization  of  Metallic  Materials.  36 

4.  Characterization  Tests  for  B/Al  and  Stainless  Steel  Wire  Mesh.  41 

5.  Parametric  Matrix  Defining  Structural  Elements  and  Impact  48 

Conditions  to  be  Analyzed. 

6.  Shape,  Size  and  Configuration  Details  of  Structural  Elements.  52 

7.  Task  VI  Fabricated  Composite  Specimen  Test  Conditions.  56 

8.  Task  VI  -  Real  Composite  Blade  Test  Conditions.  58 

9.  Task  VI  -  Fabricated  Metal  Specimen  Test  Conditions.  59 

10.  Task  VI  -  Real  Metal  Blade  Test  Conditions.  60 

11.  Matrix  of  Error  Band  Properties  and  Associated  Damage  Modes.  62 


viii 


1.0  INTRODUCTION 


The  objective  of  this  progam  is  to  develop  and  validate  structural  design 
criteria  that  account  for  the  transient  overloads  produced  by  bird  and  ice 
impacts  on  turbine  engine  first-stage  fan/compressor  blades.  This  program  is 
part  of  a  continuing  effort  of  the  Air  Force  to  minimize  ownership  costs  by 
placing  added  emphasis  on  the  derivation  of  durable  damage-tolerant  advanced 
turbomachinery  component  designs. 

Foreign  object  damage  (FOD)  in  aircraft  engines  has  been  an  increasing 
drain  on  defense  economic  resources.  Present  generation  engine  blading  de¬ 
signs  have  substantial  material  toughness  and  cross  section  to  meet  the  mechan¬ 
ical  and  aeromechanic al  requirements,  but  they  were  not  specifically  designed 
for  damage  tolerance.  The  prevailing  design  state-of-the-art  concerning 
damage-tolerant  blading  is  mostly  empirical  and  based  on  experience  from  major 
incidents  and  associated  ad  hoc  testing.  The  empirical  approach  to  designing 
damage- tolerant  blading  requires  considerable  investment  of  resources  and  is 
sometimes  incapable  of  providing  directly  predictable  foreign  object  impact 
response  of  blading.  This  approach  is  not  adequate  for  newer,  more  damage- 
prone  blading  designs  incorporating  lightweight  materials,  including  advanced 
composites,  and  thinner  cross  sections  to  achieve  improved  performance. 

To  achieve  more  efficient  damage-tolerant  blading,  comprehensive  foreign 
object  impact  design  criteria,  based  upon  transient  structural  response  tools, 
are  needed.  These  tools  will  provide  direct  assessment  of  a  blade's  impact 
damage  tolerance  and  identify  areas  for  improvement.  The  purpose  of  this  pro¬ 
gram  is  to  provide  the  necessary  computer  tools  and  validation  testing  to 
establish  reliable  foreign  object  impact  design  criteria. 

This  program  consists  of  11  tasks  which  progressively  develop  the  FOD 
design  criteria  from  computer  models,  structural  element  and  material  property 
tests,  and  static  and  rotating  single-blade  tests  to  full-  or  partial-stage 


rotating  tests. 

The  11  tasks  are  identified  as  follows: 

Task  I 

- 

Design  System  Structure 

Task  11 

- 

Transient  Response  Analysis  Models 

Task  III 

- 

Impact  Loading  Model 

Task  IV 

- 

Material  Response  and  Failure  Criteria 

Task  V 

- 

Parametric  Analysis 

Task  VI 

- 

Structural  Element  Tests 

Task  VII 

- 

Error  Band  Analysis 

Task  VIII 

- 

Foreign  Object  Impact  Design  Criteria 

Task  IX 

- 

Single-Blade  Impact  Tests 

Task  X 

- 

Full-Stage  Response  Prediction 

Task  XI 

- 

Full-Scale  Impact  Tests 

The  design  analysis  methods  and  failure  criteria  derived  in  the  course 
of  this  program  will  be  applicable  to  both  advanced  composite  materials  and 
monolithic  materials  of  construction  for  current  and  advanced  fan/compressor 
blading. 
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2.0  SUMMARY 


The  program  consists  of  11  tasks  to  establish  reliable  foreign  object 
impact  design  criteria.  The  emphasis  will  be  on  determining  the  effects  of 
the  ingestion  of  birds  and  ice  during  takeoff  and  landing,  which  present  the 
greatest  hazards  to  aircraft  engines. 

The  11  tasks  progressively  develop  the  FOD  design  criteria  from  tran¬ 
sient  structural  response  computer  models,  test  specimens  and  single-blade 
tests  to  full-  or  partial-stage  rotating  rig  tests.  The  transient  response 
computer  models  will  be  based  on  existing  mathematical  approaches  and  on  exis¬ 
ting  computing  programs  that  will  be  adjusted  as  required  for  the  specific 
application  of  this  study. 

Task  I  which  develops  a  design  structure  for  the  purposes  of  providing  a 
consistent  framework  in  trfiich  to  constitute  the  design  criteria,  guiding  the 
overall  technical  effort  and  formulating  a  data  foundation  to  define  the 
problem  environment  has  been  completed. 

Tasks  II  through  VIII  are  currently  underway  and  address  to  the  develop¬ 
ment  of  the  transient  response  and  loading  models,  local  and  structural  impact 
testing,  material  property  determination,  error  band  analysis  and  the  formu¬ 
lation  of  design  criteria. 

Initial  versions  of  the  preliminary  and  final  design  transient  structural 
response  computer  programs  have  been  developed  and  checked  out.  In  addition, 
the  coding  for  interfacing  the  bird  loading  model  developed  by  the  University 
of  Dayton  Research  Institute  with  the  final  design  response  model  developed 
by  General  Electric  has  been  completed  and  is  in  the  final  checkout  stage. 

Three  first  stage  blades  were  selected  for  modeling  and  structural  test¬ 
ing.  Composite  and  metallic  specimens  for  local  impact  and  structural  element 
bench  impact  and  material  property  tests  were  designed  and  most  of  these  speci¬ 
mens  have  been  fabricated.  Test  plans  defining  the  testing  of  these  specimens 
and  real  blades  have  been  completed  and  approved  by  the  Air  Force.  Real  blade 
impact  testing  to  establish  the  maximum  strain  rates  that  are  needed  for  the 
material  property  tests  has  been  completed. 

A  parametric  matrix  to  define  the  conditions  and  geometries  that  will  be 
analyzed  in  Tasks  V  and  X,  and  to  provide  guidance  in  the  selection  of  the 
structural  tests  that  will  be  conducted  in  Tasks  VI  and  IX,  has  been  formu¬ 
lated  . 

Work  was  started  on  the  compilation  of  relevant  blade  and  impactor  prop¬ 
erties  for  use  in  the  error  band  analysis.  FOD  design  criteria  goals  were 
established  within  the  framework  of  the  projected  capabilities  of  the  tran¬ 
sient  structural  analysis  computer  tools  being  developed  in  Tasks  II  and  III. 


3.0  TASK  I  -  DESIGN  SYSTEM  STRUCTURE 


A  design  system  structure  was  developed  through  the  Task  I  effort  to 
establish  a  consistent  framework  in  which  to  constitute  FOD  design  criteria 
and  guide  the  overall  technical  effort.  This  design  system  structure  defines 
the  objective  and  interfaces  of  the  various  multidisciplined  tasks  that  must 
be  satisfied  to  achieve  meaningful  design  criteria  and  includes  a  data  founda¬ 
tion  that  describes  the  problem  envirotment. 

Figure  l  depicts  the  interrelations  between  the  major  distinct  activi¬ 
ties  of  quantification  of  experience,  analytic  development,  bench  tests, 
rotating  rig  tests,  and  Che  work  flow  of  the  program. 

The  results  of  the  Task  I  effort  were  documented  in  a  Research/Develop¬ 
ment  Test  Plan1  that  was  approved  by  the  Air  Force. 

To  conduct  the  program  described  in  Figure  1,  a  team  of  technical  per¬ 
sonnel  has  been  gathered  from  several  organizations  to  utilize  their  specific 
technical  expertise  and  foreign  object  impact  experience.  These  organizations 
are : 

•  General  Electric  Aircraft  Engine  Group  (AEG) 

•  General  Electic  Corporate  Research  and  Development  Center  (CR&DC) 

•  University  of  Dayton  Research  Institute  (UDRI) 

Figure  2  shows  the  functional  responsibilities  of  the  participants. 


*Superscript  numbers  refer  to  references. 


Figure  1.  Program/Tasks  Flow  -  Foreign  Object  Impact  Design  Criteria. 


4.0  TASK  II  -  TRANSIENT  RESPONSE  ANALYSIS  MODELS 
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The  first-  and  second-level  models  selected  for  use  in  the  transient 
response  analysis  of  first  stage  fan  and  compressor  blades  were  approved  by 
the  Air  Force. 


4.1  INTRODUCTION 


The  transient  response  of  first  stage  fan  and  compressor  blades  due  to 
foreign  object  impact  will  be  analyzed  on  two  different  levels  during  the 
course  of  this  program.  The  rationale  for  this  approach  is  based  on  the  fact 
chat  each  of  the  levels  of  analysis  possess  intrinsic  advantages  which  make 
them  very  attractive  for  the  study  of  different  aspects  of  the  FOD  problem. 

The  first  level  of  analysis  will  be  accomplished  with  a  modified  version  of 
the  finite  element  program  NONSAP.  This  analysis  will  be  capable  of  carrying 
out  a  very  detailed  analysis  embracing  local  impact  damage  as  well  as  gross 
structural  damage.  NONSAP  is  a  very  general  and  accurate  package  for  non¬ 
linear,  finite  element  analyses  and  contains  3-D  elements  trfiich  will  be  used 
herein.  The  modified  version  which  will  be  used  in  this  program  will  be  cal¬ 
led  NOSAPM.  This  program  has  a  library  of  three  different  time  integration 
techniques  (Wilson  e,  Newmark,  central  difference)  as  well  as  two  elasto- 
plastic  hardening  models  (kinematic  hardening  and  isotropic  hardening).  This 
versatility  will  prove  useful  in  finding  the  most  efficient  strategy  for  solu¬ 
tion  of  the  dynamic  problem.  However,  the  detail  and  versatility  of  this  NON¬ 
SAP  package  is  obtained  at  the  expense  of  high  computing  costs  tfiich  could 
prove  to  be  a  problem  during  parametric  design  studies.  As  a  result,  the 
second  level  analysis  based  on  S.  Levy's  component  element  method  of  tran¬ 
sient  dynamics^  will  be  employed  in  these  situations  where  cost  becomes  an 
important  factor.  The  savings  in  the  Component  Element  Method  (COMET) 
are  accomplished  through  the  use  of  simplified  element  models  and  an 
explicit  time  integration  technique  which  allows  nonlinear  structural  dynamics 
problems  to  be  solved  without  ever  physically  assembling  stiffness  matrices. 
This  Level  2  analysis  will  be  somewhat  less  detailed,  but  experience  in 
other  fields  including  pipe  whip4,  locomotive  dynamics,  and  nuclear  fuel 
rod  vibration^  shows  that  the  technique  is  an  accurate  as  well  as  efficient 
tool.  In  addition  to  being  of  use  in  parametric  studies,  the  lower  cost  of 
this  level  of  analysis  should  prove  to  be  an  advantage  in  blade-interaction 
studies . 

No  dynamic  crack  growth  analyses  will  be  incorporated  into  the  first  and 
second  level  analyses.  Stresses  and  strains  will  simply  be  calculated  and 
then  used  to  indicate  damage  based  on  the  design  criteria  to  be  formed  as  part 
of  Task  VIII.  Plastic  curlback  of  blades  will  be  considered.  Graphical  output 
of  stresses  and  strains  will  be  provided,  and  a  method  of  implementing  the 
results  of  the  impact  damage  design  criteria  within  the  structural  codes  will 
be  incorporated.  Based  on  the  transient  response  models  of  Task  II  and  the 
damage  criteria  formulated  in  Task  VIII,  the  final  package  would  then  indicate 
probable  areas  of  failure  as  they  are  encountered  in  the  analysis. 
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4.2  FIRST-LEVEL  RESPONSE  MODEL  (NOSAPM) 


As  indicated  above,  the  program  NOSAPM  contains  a  vide  range  of  tools 
vhich  will  be  necessary  in  modeling  the  FOD  problem.  Specifically,  in  this 
analysis,  the  3-D  isoparametric  element  will  be  used.  This  element  can  be  used 
with  a  variable  number  of  nodal  points  ranging  from  8  to  21.  It  is  antici¬ 
pated  that  either  16  noded  or  20  noded  elements  will  be  used  to  allow  the 
blades  to  be  modeled  with  one  element  through  the  thickness  tfcile  still  in¬ 
cluding  the  bending  effects  which  will  be  most  important  in  this  problem. 

NOSAPM  is  capable  of  the  large-displacement,  elasto-plastic  analysis,  and  uses 
an  out-of-core  equation  solver  for  efficient  solution  of  large  problems  such 
as  this  one.  The  inelastic  material  behavior  will  be  modeled  with  a  bilinear 
stress-strain  curve  and  both  isotropic  and  kinematic  hardening  will  be  avail¬ 
able  for  application.  In  the  case  of  composite  materials,  a  linear  ortho¬ 
tropic  material  model  is  already  available.  If  metal-matrix  composites  must 
be  analyzed  and  anisotropic  yielding  is  of  concern,  it  appears  that  this  effect 
can  be  accounted  for  with  the  use  of  Hill's  theory  of  anisotropic  yielding. 
Although  not  presently  programmed  in  NOSAPM,  the  level  of  experience  here  in¬ 
dicates  that  it  would  not  be  overly  difficult  to  accomplish.  Centrifugal 
stiffening  has  been  incorporated  in  NOSAPM  through  the  addition  of  a  body- 
force  (centrifugal)  load  vector.  This  load  vector  is  used  to  calculate  cen¬ 
trifugal  displacements  which  are  used  as  initial  conditions  for  the  transient 
response.  The  constant  centrifugal  load  vector  is  also  added  to  the  impact 
loads  for  each  time  step.  The  centrifugal  displacements  calculated  by  the 
initial  state  of  the  blade  before  the  impact  are  also  utilized  in  the  calcu¬ 
lation  of  eigenfrequencies  and  eigenmodes  that  indicate  the  effects  of 
centrifugal  stiffening. 

If  strain-rate  effects  are  shown  to  be  substantial  in  the  course  of 
material  characterization,  they  will  also  be  included.  Flexible  root  fixity 
as  well  as  blade  interaction  at  the  shrouds  will  be  handled  with  spring 
elements.  Figures  3(a)  and  (b)  show  the  finite  elements  of  the  blade. 

A  successful  compilation  of  the  current  version  of  the  NOSAPM  source 
file  placed  on  the  WPAFB  CDC6600  computer  system  has  been  accomplished 
through  the  remote  I/O  site  at  Evendale.  In  addition,  a  successful  check¬ 
out  problem  was  run.  This  consisted  of  the  fiat  plate  shown  in  Figure  4, 
modeled  with  20  noded  elements,  rotating  at  5000  rpm  and  loaded  with  a 
transient  load . 

In  order  to  facilitate  the  running  of  NOSAPM  from  Evendale,  a  data  deck 
generator  program  is  being  written.  This  program  is  90Z  complete  and  will 
generate  card  images  of  the  NOSAPM  input  data  based  on  element  and  node  files 
(created  with  an  existing  mesh  generator  program  available  at  Evendale)  and 
on  user  input  specified  in  the  interactive  time-sharing  mode.  The  input 
data  decks  so  created  will  be  transmitted  to  the  WPAFB  6600  systems  via  the 
remote  Evendale  1/0  site  to  establish  data  files.  Control  of  job  execution 
in  the  batch  mode  is  then  accomplished  at  Evendale  through  time-sharing  with 
the  job  output  directed  to  the  remote  I/O  site  at  Evendale. 


Three-Dimensional 
Finite  Elements 


Level  1  Response  Model  Formulated 
With  3-D  Finite  Elements. 
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4.3  SECOND-LEVEL  RESPONSE  MODEL  (COMET-BLADE) 


Essentially,  the  same  physical  effects  are  considered  in  the  Level  2 
(COMET -BLADE)  analysis,  but  the  elemental  approach  is  somewhat  simplified, 
thus  contributing  to  gains  in  efficiency.  The  Level  2  analysis  is  based 
on  the  component  element  method^  and  provides  an  analytical  computer  capa¬ 
bility  for  predicting  gross  and  local  structural  response  of  blades  to  re¬ 
presentative  large  body  impact  loading  stimuli. 

Bending  and  membrane  effects  are  modeled  with  a  system  of  plate  elements. 
These  elements  couple  membrane  and  bending  forces  and  are  similar  to  those 
that  have  been  highly  successful  in  modeling  both  plasticity  and  hysteresis 
in  the  General  Electric  COMET-PIPE  computer  program^. 

The  formulation  and  programming  of  the  plate  elements  for  large  deforma¬ 
tions,  material  yielding  and  hysteresis  effects  are  complete.  These  elements 
have  nine  nodes  and  address  to  three  cases:  (a)  an  element  around  a  central 
node  for  interior  regions  of  the  blade;  (b)  a  corner  element;  and  (c)  a  raiddle- 
of-an-edge  element.  Figure  5  gives  an  indication  of  the  Level  2  modeling. 

The  Level  2  analysis  also  includes  centrifugal  loading  effects  in  a 
manner  similar  to  that  of  the  first  level  analysis  and  provides  for  the  speci¬ 
fication  of  damping  that  can  be  proportional  to  mass,  to  stiffness,  or  to  a 
combination  of  the  two.  The  blade  can  have  twist  and  camber  and  can  vary  in 
thickness.  Boundary  conditions,  such  as  restraint  of  rotation  and  displace¬ 
ment,  can  be  specified. 

The  blade  is  modeled  by  giving  the  positions  of  nodes  in  global  coordi¬ 
nates,  the  blade  thickness,  the  material  real-stress  versus  real-strain  curve, 
the  density,  the  damping,  the  speed  in  rpm,  and  the  boundary  conditions. 

The  loads  can  be  applied  at  any  node  in  any  direction.  Step-by-step  numerical 
integration  is  used  to  obtain  a  solution  for  the  response.  Still  to  be  in¬ 
corporated  into  the  second  level  analysis  is  a  material  model  that  addresses  to 
linear  anistropic  effects  and  spring  elements  that  allow  for  blade-to-blade 
coupling.  This  latter  feature  should  be  very  useful  in  modeling  the  inter¬ 
action  of  several  blades  due  to  the  advantage  of  lower  computing  costs  for 
the  COMET  analysis. 


4.4  INTERFACE  WITH  THE  LOADING  MODEL 


The  interface  between  the  response  and  loading  models  will  be  described 
in  detail  in  the  loading  model  section  (Section  5) .  A  brief  description  of  the 
interface  follows. 


Since  the  impact  loading  due  to  the  foreign  object  will  depend  upon  the 
deformation  of  the  blade,  interaction  between  the  structural  model  and  the 
impactor  model  will  be  necessary.  This  requires  the  interfacing  of  the 
response  routines  developed  by  CR&D  with  the  loading  model  developed  by  UDRI. 
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The  interface  between  the  first-level  response  model  developed  by  CB&D  and 
the  loading  model  developed  by  UDR1  has  been  accomplished.  Geometric  and 
dynamic  information  is  periodically  transferred  from  the  structural  analysis 
and  used  to  update  the  impact  loads.  This  need  not  be  done  at  every  structural 
time  step  of  the  transient  analysis  since  large  time  steps  may  be  suitable 
without  loss  of  accuracy.  The  interface  with  the  loading  model  will  also  be 
included  in  the  second  level  response  model.  Figure  6  shows  the  linkup 
between  the  loading  model  subroutine  and  the  executive  or  main  response  model 
program. 
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Figure  6.  Link-up  Between  the  Loading  Model  Subroutine 

and  the  Executive  or  Main  Response  Model  Program. 
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5.0  TASK  III  -  IMPACT  LOADING  MODEL 


5 . 1  INTRODUCTION 


The  objeccive  of  chis  cask  is  Co  geaerace  an  impact  loading  model  for 
birds  and  ice  which  can  be  inCerfaced  and  used  wich  Che  Cransienc  response 
sCrucCural  analysis  Cools  developed  in  Task  II.  The  efforc  in  chis  Cask  is 
divided  inco  four  subcasks:  (a)  sioulaced  bird  and  ice  impacc  parameters; 
(b)  loading  model  sensitivity  study;  (c)  formulae ion/ incer facing  of  loading 
model  with  structural  models;  and  (d)  experimental  verification  of  loading 
mode  1 . 


5 . 2  SUBTASK  A  -  SIMULATED  BIRD  AND  ICE  IMPACT  PARAMETERS 


The  work  on  chis  subcask  has  been  completed.  A  report  written  by  UDRI 
documenting  the  effort  is  currently  being  reviewed  preparatory  to  trans¬ 
mittal  to  the  Air  Force.  This  report  presents  the  results  obtained  from 
spatially  and  temporally  resolved  pressure  measurements  made  during  impacts 
of  gelatin  artificial  birds  and  ice  on  rigid  targets.  Data  are  presented  that 
show  that  gelatin  behaves  the  same  as  real  birds  during  impact  and  that  the 
analytically  predicted  pressure  distributions  obtained  with  the  loading  model 
of  Subtask  C  agree  with  the  measured  pressures  for  both  bird  and  ice 
impacts . 


5 . 3  SUBTASK  B  -  LOADING  MODEL  SENSITIVITY  STUDY 


The  purpose  of  the  sensitivity  study  is  to  provide  specific  guidance 
as  to  which  aspects  of  the  loading  model  formulation  must  be  modeled  most 
accurately  to  properly  predict  specific  aspects  of  blade  response.  The  MARC/CDC 
finite-element  code  is  being  used  for  the  numerical  analysis  in  the  sensitivity 
study.  Nine,  twenty-four  and  thirty-two  20-node  brick  element  models  were 
used  to  investigate  the  sensitivity  of  the  response  of  typical  blade-like 
structures  to  the  details  of  impact  loading.  The  9-  and  24-element  models 
were  used  to  predict  the  response  of  a  flat  cantilevered  place  with  the 
aspect  ratio  of  the  F101  first  stage  fan  blade;  Figure  7  shows  the  9-element 
model.  The  results  obtained  wich  these  9-  and  24-element  models  indicate  chat 
the  momentum  transferred  to  the  plate  is  of  primary  importance.  No  signifi¬ 
cant  local  deformations  were  observed  during  loading  with  these  models.  The 
32-element  model  was  used  to  determine  the  importance  of  spatial  and  temporal 
loading  on  the  local  deformations  of  plate-like  structures  with  nonconstant 
cross  section.  This  model  provided  the  relationship  between  loading  (includ¬ 
ing  the  importance  of  initial  high  peak  pressures)  geometry,  and  mesh  size  on 
a  detailed  basis. 
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All  of  the  required  computer  runs  have  been  completed  and  the  results 
have  been  factored  into  the  formulation  used  for  the  loading  model.  A  report 
summarizing  the  results  is  nearing  completion. 


5.4  SUBTASK  C  -  FORMULATION/ INTERFACING  OF  LOADING  MODEL  WITH  STRUCTURAL  MODELS 

5.4.1  Introduction 

The  loads  exerted  on  a  structure  during  an  impact  can  be  characterized  in 
terms  of  a  number  of  fundamental  parameters.  The  most  important  of  these 
parameters  is  the  total  impulse  imparted  to  the  structure,  the  forces  (i.e., 
the  rate  of  impulse  transfer),  the  pressures  (i.e.,  the  spatial  distribution 
of  the  force),  and  the  time  variation  of  forces  and  pressures.  The  parameters 
that  must  be  considered  in  a  particular  impact  are  a  function  of  the  response 
regimes  of  the  target.  We  initially  identified  three  response  regimes  and 
described  suitable  loading  models  for  each  regime.  The  Level  1  loading  model 
was  applicable  to  rigid,  immovable  targets.  The  Level  2  loading  model  was 
applicable  to  targets  that  translate  and  rotate  but  do  not  deform  locally. 

The  Level  3  loading  model,  the  most  general  of  the  three,  was  applicable  to 
impacts  in  which  the  target  not  only  translates  and  rotates  during  impact, 
but  also  deforms  locally. ^ 

During  the  course  of  the  present  contract,  we  decided  that  only  the 

Level  3  loading  model  was  of  major  interest  because  of  the  response  regime  of 

fan  blades  and  the  relative  ease  of  implementation.  Subsequently,  the  orig¬ 
inal  Level  3  loading  model  was  renamed  as  the  Level  1  loading  model.  Level 
1  was  thus  used  to  denote  both  the  most  general  loading  model  and  the  most 
general  structural  analysis  model. 

Described  herein  is  the  development  of  the  Level  1  loading  model  for  pre¬ 
dicting  the  pressure  distribution  on  a  surface  impacted  by  a  bird  or  ice.  The 

model  is  based  on  a  fluid  dynamic  description  of  soft  body  impacts.  Experi¬ 

mental  evidence  in  tue  form  of  measured  pressure  distributions  (spatial  and 
temporal)  during  impact  strongly  supports  the  concept  that  bird  and  ice  impacts 
are  basically  fluid  dynamic  in  nature.^ 

The  loading  model  developed  specifically  addresses  the  problem  of  birds, 
ice  spheres,  and  slab  ice  impacting  fan  and  compressor  blades.  The  overall 
model  consists  of  a  model  for  the  slicing  of  the  impacting  material  by  a  row 
of  blades  and  of  a  fluid  dynamic  model  for  describing  the  impact  of  an  indivi¬ 
dual  slice  on  a  blade.  The  model  is  set  up  to  interface  with  the  structural 
analysis  computer  programs  developed  in  Task  II. 

Described  below  are  the  slicing  model,  Che  fluid  flow  model  of  impact, 
the  structure  of  the  loading  model  subroutine,  and  the  method  used  to  inter¬ 
face  the  loading  model  subroutine  with  the  structural  response  computer  pro¬ 
grams  . 
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5.4.2  The  Slicing  Model 

Because  Che  forces  required  Co  slice  a  sofc  body  are  small  compared  Co 
Che  forces  required  Co  decelerace  Che  body,  Che  slicing  forces  can  be  ignored. 
Therefore,  Che  slicing  model  developmenC  reduces  Co  simply  a  geometric  problem 
of  deCermining  Che  dimensions  and  veighc  of  a  slice.  For  birds  and  ice 
spheres  (such  as  hail),  Che  velocity  of  Che  ingesCed  objecC  is  small  compared 
Co  the  velociCy  of  Che  aircraft;  thus,  che  object  velocity  can  be  ignored. 

For  slab  ice,  such  as  ice  breaking  loose  from  an  engine  nacelle,  Che  velocity 
of  che  slab  relacive  Co  che  nacelle  is  not  well  defined;  buc  we  can  expecc 
chat  velociCy  is  Co  be  small  in  comparison  Co  Che  aircrafc  velocity. 


5. 4. 2.1  Slicing  Model  Development  for  Birds 

The  bird  is  idealized  as  a  righc  circular  cylinder  with  a  length-Co- 
diameter  ratio  of  2,  and  che  velocity  of  Che  bird  relacive  Co  che  aircrafc 
is'caken  to  be  equal  and  opposite  Co  Che  aircrafc  velocity. 

The  bird/blade  interaction  geometry  is  shown  in  Figure  8.  The  given 
inpuC  daCa  muse  include  Che  number  of  blades  per  stage,  Che  rotational  speed 
of  Che  blade  (rpm),  Che  radius  from  Che  rotational  axis  of  Che  rotor  Co  che 
point  on  che  blade  at  which  Che  center  of  Che  impact  occurs,  and  che  blade 
orientation  angle. 

An  infinite  number  of  slice  shapes  is  possible  for  a  given  set  of  these 
inpuC  parameters;  che  shape  depends  on  Che  orientation  of  Che  bird  relative 
Co  Che  blade  and  on  Che  span  location  on  Che  blade  at  which  che  impacC  occurs. 
Since  worse-case  slice  shapes  are  desired  (i.e.,  slice  shapes  having  Che 
largest  possible  slice  mass),  Che  orientation  of  the  axis  of  Che  bird  and  the 
center  of  impacc  of  a  slice  were  chosen  to  produce  a  slice  having  Che  largest 
possible  mass.  Figure  8  illustrates  this  orientation  for  a  bird  idealized  a; 
a  right  circular  cylinder  with  a  length-to-diameter  ratio  of  2.  The  largest 
slice  mass  occurs  when  Che  axis  of  che  slice  is  coincidenC  with  Che  axis  of 
che  righc  circular  cylinder.  The  bird  slice  mass  is  given  by  che  following 
equation: 

Bird  Slice  WeighC  «  W8  ■ 

-  20D  J\^(D/2)2  -  (h/2)2+  ~  sin-1 

Where  p  »  bird  density  *  0.91  g/cm^ (0.033047  lbf/in^)  and  the  other  terms  are 
defined  in  Tables  1  and/or  2. 

The  bird  weight  and  Che  bird  density  together  with  che  slice  shape  define 
che  slice  mass.  The  amount  of  bird  slice  consumed  ac  any  time  during  impact  is 
directly  related  Co  Che  velocity  of  Che  bird  relacive  Co  che  blade  (Vr  in 
Figure  8)  and  Che  Cime  from  che  beginning  of  impacc. 
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-Ulatle/Interaction  Geometry. 


Blade  Calculation  Input  Data  Sheet  lor  3  Oz  Bird 
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Table  2.  Blade  Calculation  Input  Data  Sheet  for  1.5-lb  Bird. 
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In  addition  Co  the  idealized  right  circular  cylinder  slicing  model  for 
birds,  slicing  models  will  be  also  developed  for  ice  balls  and  slab  ice. 


5. 4. 2. 2  Slice  Geometry  Parametric  Study 

Bird  slice  parameters  based  on  the  right  circular  cylinder  bird  have 
been  computed  for  both  starling  impacts  ( 3-ounce  weight)  and  larger  bird 
impacts  (1.5-pound  weight).  Three  different  blade  configuration  types 
(J79,  APSI,  and  F101  engines)  were  considered;  and  typical  values  of  the 
rotor  speed,  number  of  blades  per  stage,  blade  orientation  angle,  and  impact 
radius  were  used. 

Tables  1  and  2  show  the  input  and  computed  parameters  for  the  two  bird 
size  impacts.  The  primary  observation  from  this  parameter  study  is  that 
these  impacts  are  highly  oblique  (impact  angles  ranging  from  51°  down  to  19°). 
The  bird  slice  model  has  been  incorporated  into  the  loading  model.  For  the 
loading  model,  the  bird  slice  model  essentially  defines  the  velocity  and 
shape  of  a  fluid  jet  impacting  a  blade  as  well  as  the  duration  of  the  impact. 


5.4.3  Fluid  Flow  Model  of  Impact 

It  has  been  demonstrated  in  Reference  6  that  for  birds  having  an  L/D  of  2 
or  greater  the  largest  contributing  factor  in  the  total  impulse  of  a  bird 
impacting  a  target  is  the  steady  flow  phase  of  an  impact.  Thus,  in  the  cur¬ 
rent  version  of  the  loading  model,  the  pressure  distribution  is  based  on  a 
steady  flow  analysis.  This  method  of  analysis  is  extended  into  a  dynamic 
analysis  for  deforming  targets  by  treating  the  flow  at  any  instant  of  time  as 
a  quasi-steady  flow  using  the  instantaneous  relative  velocity  between  the  bird 
slice  and  the  deforming  target  as  the  characteristic  velocity.  The  instan¬ 
taneous  shape  of  target  is  used  to  define  the  surface  on  which  impact  occurs - 

In  the  analysis  of  the  impact  process,  the  strength  of  the  projectile 
material  is  considered  to  be  negligible.  This  assumption  is  quite  reasonable 
for  the  typical  projectile  and  target  materials  of  interest. 

We  assume  that  at  least  as  a  first  approximation,  the  fluid  flow  may  be 
treated  as  incompressible.  This  is  a  reasonable  assumption  because  the 
measured  steady-state  pressures  are  quite  small  in  comparison  to  the  pressures 
required  to  produce  significant  density  changes  in  birds  or  ice. 

By  virtue  of  the  incompressible  flow  and  steady  flow  assumptions  and  the 
negligible  strength  of  the  projectile  material,  the  problem  of  predicting  soft 
body  impact  loads  is  amendable  to  analysis  within  the  following  real  con¬ 
straints.  The  computing  time  required  to  obtain  a  reasonable  solution  must  be 
small  in  comparison  to  the  computing  time  required  to  compute  the  structural 
response,  and  the  computer  storage  requirement  for  the  loading  model  must  be 
small  in  comparison  to  the  structural  response  analysis  computer  program 
storage  requirements. 
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Three-dimensional  potential  flow  theory  was  chosen  as  the  most  appro¬ 
priate  method  for  modeling  the  impact  process.  Our  initial  work  in  using  this 
approach  was  quite  successful.  In  that  work,  we  used  the  method  of  surface 
singularities  to  determine  the  velocity  and  pressure  fields  caused  by  circular 
jets  impacting  at  oblique  angles  on  flat  rigid  plates.  This  same  technique 
was  used  to  develop  the  loading  model.  The  method  developed  earlier  was 
generalized  to  include  the  impact  surface  curvature  and  deformation  velocity, 
an  arbitrary  cross-sectional  area  of  the  impacting  flow,  and  a  method  for 
generating  planar  quadrilateral  surface  elements  from  given  finite-element 
surface  nodal  point  locations. 


5. 4. 3.1  Surface  Singularity  Method  -  Superposition  of  Onset 
Flow  and  Surface  Source  Distribution 


r 


The  impacting  bird  slice  is  idealized  as  an  onset  flow  of  velocity,  Vr, 
uniformly  distributed  over  the  cross-sectional  area  of  the  slice.  The  fluid 
velocity  at  a  point  is  expressed  as  the  negative  gradient  of  a  potential 
function  that  satisfies  Laplace’s  Equation  in  the  region  exterior  to  the 
blade,  that  has  a  zero  normal  derivative  on  the  blade  surface,  and  that 
approaches  the  proper  uniform  stream  potential  function.  The  flow  is  viewed 
as  an  onset  flow  potential  and  disturbance  potential  due  to  the  blade.  The 
disturbance  potential  satisfies  the  Laplace’s  Equation  in  the  exterior  and 
the  appropriate  normal  derivative  condition  on  the  surface  and  vanishes  at 
infinity.  From  potential  theory,  it  can  be  shown  that  this  potential  may  be 
evaluated  in  terms  of  a  surface  source  density  distribution  with  which  the 
blade  surface  may  be  considered  to  be  covered.  The  incorporation  of  the 
normal  derivative  condition  on  the  surface  produces  a  source  density  dis¬ 
tribution  described  by  a  two-dimensional  Fredholm  Integral  Equation  of  the 
second  kind  over  the  surface.  Once  this  equation  is  solved  for  the  source 
density  distribution,  first  the  disturbance  potential  and  then  the  distur¬ 
bance  flow  velocities  may  be  obtained. 


5. 4. 3. 2  The  Body  Surface  Approximation 

In  the  numerical  solution  of  the  Fredholm  Integral  Equation  for  the  sur¬ 
face  source  density  distribution,  the  body  surface  must  be  suitably  approxi¬ 
mated.  The  body  surface  is  approximated  by  a  large  number  of  small  plane 
elements  that  are  formed  from  the  original  points  defining  the  body  surface. 
The  calculations  employ  plane  quadrilateral  elements,  and  the  source  density 
is  assumed  constant  over  each  of  these  elements.  With  this  assumption,  the 
problem  of  determining  the  continuous  function  for  the  source  density  reduces 
to  one  of  determining  a  finite  number  of  values  -  one  for  each  of  the  planar 
elements . 

5. 4. 3. 3  Solution  Procedure 

This  section  presents  a  brief  description  of  the  method  of  solution. 
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As  a  first  step,  it  is  necessary  to  prescribe  the  input  points,  a  set 
of  points  in  three-dimensional  space  describing  the  body  surface.  These 
input  points  should  be  chosen  to  provide  the  best  representation  of  the 
body  with  the  fewest  possible  points.  Four  of  these  inputs  are  employed  to 
generate  one  plane  quadrilateral  source  element  of  the  body  surface. 

Once  the  values  of  the  surface  source  density  have  been  determined,  the 
velocity  components  at  each  null  point  are  calculated;  and  the  appropriate 
components  of  the  onset  flow  are  added.  The  resultant  velocities  and  pres¬ 
sures  are  easily  computed  from  the  velocity  components. 

5. 4. 3. 4  Onset  Flow  Modeling 

When  it  is  used  with  the  surface  singularity  technique,  the  simple  onset 
flow  model  described  in  Section  5. 4. 3.1  is  adequate  for  normal  impacts  and 
also  for  oblique  impacts,  provided  that  the  angle  of  incidence  (the  angle 
9  in  Figure  8)  is  greater  than  about  30°.  For  small  values  of  9,  a  more 
sophisticated  onset  flow  model  is  used.  The  basis  for  this  more  accurate 
onset  flow  model  is  a  two-dimensional  jet  flow  theory  for  oblique  impacts. 


5.4.4  Structure  of  the  Loading  Model  and  Interface  With  Response  Models 

We  have  incorporated  in  our  analysis  a  procedure  for  finite-difference 
numerical  computation  of  the  potential  flow  about  arbitrary  three-dimensional 
bodies  that  has  certain  features  that  render  it  compatible  with  the  finite- 
element  structural  response  analysis  code.  Specifically,  the  nodal  points  in 
the  finite-element  program  become  approximate  locations  of  the  null  point 
(defined  to  be  the  point  where  the  quadrilateral  element  induces  no  velocity 
in  its  own  plane)  around  which  the  plane  quarilateral  elements  are  constructed 
in  the  potential  flow  analysis.  This  greatly  facilitates  the  coupling  of  the 
finite-difference  loading  model  and  the  finite-element  structural  code. 

In  loading  used  for  any  arbitrarily  shaped  impact  surface,  the  impact 
area  is  divided  into  small  flat  elements;  and  a  uniform  distribution  of 
sources  is  assumed  to  cover  each  area.  At  the  beginning  of  impact,  an  initial 
pressure  distribution  is  computed  for  the  undeformed  blade.  During  an  impact 
in  which  local  deformation  takes  place,  the  deformed  shape  of  the  impact  zone 
is  calculated  in  the  dynamic  structural  response  analysis.  After  significant 
deformation  has  occurred,  the  geometry  of  the  impact  zone  is  provided  to  the 
loading  model.  The  loading  model  is  then  used  to  calculate  a  new  pressure 
distribution.  As  the  structural  analysis  calculation  proceeds,  the  local 
shape,  the  location,  and  the  velocity  of  the  impact  area  are  updated  and 
passed  to  the  loading  model  at  appropriate  time  intervals.  The  loading  model, 
in  turn,  provides  updated  pressure  distribution  information  for  the  struc¬ 
tural  response  computation.  The  loading  model  is  fully  interactive  with  the 
structural  response  calculation.  The  duration  of  the  impact  is  computed  by 
keeping  track  of  how  much  slice  has  been  consumed  during  each  time  interval. 


A  summary  of  the  interface  between  che  loading  and  response  models  fol¬ 
lows: 

•  The  response  model  provides  displacements  and  velocities  at  the  sur¬ 
face  nodes  of  the  structural  elements  in  the  impact  region.  These 
structural  elements  are  then  subdivided  into  smaller  elements  corre¬ 
sponding  to  the  finer  grid  of  the  loading  model  as  specified  by  user 
input,  and  the  structural  element  shape  functions  are  then  used  to  ‘ 
compute  the  displacements  and  velocities  at  the  nodes  of  the  smaller 
elements.  These  displacements  and  velocities  are  provided  to  the 
loading  model  at  specified  time  steps. 

•  The  slice  size  calculations  are  based  on  user  input  and  are  per¬ 
formed  by  the  response  program.  The  impactors  that  will  be  avail¬ 
able  to  the  user  include  the  following: 

-*  Bird  of  any  weight  represented  by  a  right  circular  cylinder 
having  L/D  ratio  equal  to  2.0. 

-  Ice  ball  of  any  diameter. 

-  Slab  ice  of  any  thickness  and  length. 

User -provided  information  (in  addition  to  the  geometry  and  mate¬ 
rial  properties  of  the  finite  element  response  model)  will  consist 
of  the  velocity  of  the  aircraft,  the  rotational  velocity  of  the  fan 
or  compressor,  the  number  of  blades  in  the  stage,  the  spanwise  loca¬ 
tion  of  the  impact,  and  the  specif ication  of  the  impactor. 

•  The  loading  model  provides  Che  pressures  at  the  nodes  of  the  smaller 
elements,  and  the  structural  element  shape  functions  are  then  used 
to  calculate  the  pressures  at  the  nodal  points  of  the  structural 
elements.  These  pressures  are  then  converted  to  equivalent  nodal 
focus  through  a  consistent  force  formation. 

•  A  flag  is  returned  by  the  loads  model  when  the  impactor  has  been 
consumed  (i.e.,  when  it  has  traveled  its  own  length).  When  this 
occurs,  the  response  model  will  no  longer  call  the  loading  model 
subroutine  and  the  nodal  pressures  are  set  to  zero. 

The  loading  model  is  capable  of  detailed  interaction  with  the  structural 
response  model  and  of  dealing  with  target  translation,  rotation,  and  local 
deformation.  The  load/response  coupling  modeled  in  this  formulation  should 
be  capable  of  accurately  predicting  both  overall  target  response  and  local 
deformation. 

The  principal  limitation  of  the  loading  model  is  that  is  does  not  include 
transient  effects.  The  most  significant  transient  effect  is  the  shock  effect; 
however,  the  porosity  present  in  birds  appreciably  reduces  the  shock  stresses 
without  significantly  affecting  steady  flow  pressures.  In  addition,  impact 
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obliquity  reduces  the  relative  importance  of  shock  stresses.  Therefore,  it  is 
not  obvious  that  neglect  of  the  shock  aspects  of  bird  impact  on  blades  is  a 
significant  deficiency.  Another  transient  aspect  of  bird  impacts  is  the  vari¬ 
ation  of  slice  cross  section  at  the  target  surface  during  impact.  This  varia¬ 
tion  results  mainly  in  a  time  variation  of  the  impact  area.  If  the  flow 
remains  quasi-steady  during  these  variations  (i.e.,  the  "velocity"  of  the 
variation  is  low  with  respect  to  local  sound  speed  -  probably  a  good  assump¬ 
tion),  then  the  loading  model  can  be  modified  to  described  these  effects.  The 
size  and  geometry  of  the  "jet"  that  flows  onto  the  target  surface  must  be  up¬ 
dated  incrementally  to  describe  the  variation  with  time  of  the  impact  area. 


5.4.5  Summary 

A  loading  model  has  been  developed  for  soft  body  impacts  on  structures 
that  respond  to  the  impact  in  a  very  general  manner.  The  loading  model  was 
developed  from  a  fluid  dynamic  description  of  a  soft  body  impacting  (flowing 
onto)  a  structure.  The  structure  can  translate,  rotate,  and  deform  locally 
during  the  impact.  The  model  is  based  on  an  approximate  solution  to  the  pro¬ 
blem  of  the  quasi-steady  flow  of  a  jet  of  arbitrary  cross  section  flowing  onto 
a  surface  with  an  arbitrary  shape. 

The  cross-sectional  area  of  the  jet  is  determined  from  an  analysis  of  the 
slicing  of  soft  bodies  by  a  row  of  fan  or  compressor  blades.  The  loading 
model  pressure  distributions  show  excellent  agreement  with  pressure  distribu¬ 
tions  obtained  in  experiments  in  which  birds,  substitute  birds,  and  ice 
spheres  impacted  on  rigid  targets.  The  loading  model  has  been  successfully 
interfaced  with  the  first  level  finite  element  structural  analysis  computer 
program. 

The  overall  validity  and  limitations  of  the  model  will  be  determined  by 
comparing  the  structural  response  predicted  by  the  combined  loading  model/ 
structural  analysis  model  with  that  observed  in  experiments. 


5.5  SUBTASK  D  -  EXPERIMENTAL  VERIFICATION  OF  LOADING  MODEL 


In  order  to  check  the  adequacy  of  the  integrated  loading  and  response 
models,  comparisons  will  be  made  between  the  analysis  predictions  and  the  test 
results  for  impacts  on  cantilevered  plates.  A  range  of  impact  parameters 
(mass  for  substitute  birds  and  ice,  velocity  and  incidence)  will  be  employed 
to  ensure  a  full  range  of  specimen  response  modes  (i.e.,  elastic,  elastic- 
plastic  structural  deformation,  and  severe  local  deformation).  The  predic¬ 
tions  of  the  complete  model  will  be  compared  to  experimentally  derived  results 
obtained  from  Task  VI  -  Structural  Element  tests.  For  some  of  cantilever 
plate  tests  of  Task  VI,  a  Moire  fringe  technique  will  be  used  to  obtain  dis¬ 
placements.  We  will,  therefore,  accomplish  the  objectives  of  both  Task  III-D 
and  VI  with  one  set  of  experiments. 
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Once  che  selected  structural  element  impact  tests  are  accomplished,  a 
detailed  comparison  between  the  experimental  and  analytic  results  will  be  made 
and  the  source  of  any  discrepancies  will  be  identified.  A  modifications 
required  in  the  loading  model  which  are  identified  in  thi'  suotask  will  be 
made. 


6.0  TASK  IV  -  MATERIAL  RESPONSE  AND  FAILURE  CRITERIA 


6.1  SUBTASK  A  -  GROSS  STRUCTURAL  DAMAGE  PROPERTIES 


6.1.1  Introduction 

Testing  is  proceeding  along  two  fronts  in  order  to  establish  gross  struc¬ 
tural  material  response  and  failure  data.  First,  bench  impact  tests  have  been 
conducted  using  the  full-scale  blades  selected  in  Task  V.  The  strain  rates 
pertinent  to  the  problem  of  foreign  object  impacts  have  been  determined  from 
these  tests.  Impact  tests ;  using  che  parameters  of  Tables  1  and  2,  have  been 
conducted  on  three  different  type  blades,  two  of  which  are  unshrouded.  The 
two  unshrouded  blades  are  the  stainless  steel  J79  first  stage  compressor  blade 
and  the  B/Al  composite  APSI  first  stage  fan  blade.  These  blades  were  canti¬ 
lever-mounted  at  the  root.  The  third  blade  type  is  Che  titanium  F101  first 
stage  fan  blade.  This  blade  is  tip  shrouded  and,  therefore,  was  fastened  at 
both  ends.  Each  blade  type  was  impacted  at  the  30  and  70%  span  locations 
and  two  blades  of  each  blade  type  were  used. 

In  the  second  part  of  Subtask  A,  tests  are  being  conducted  to  obtain 
material  property  data  for  the  relevant  strain  rates.  Forged  bar  stock  was 
used  to  make  the  metal  specimens  (8AI-IM0-IV  titanium  to  represent  the  F101  1 
blade,  and  410  stainless  steel  to  represent  the  J79  blade).  The  410  is  being 
substituted  for  the  403  which  is  unavailable.  Chemically,  the  two  alloys  are 
the  same  except  for  a  small  amount  of  chromium,  an  element  which  does  not 
affect  the  room  temperature  mechanical  properties.  Standard  cylindrical 
specimens  are  used  for  the  tension  teats  at  low  and  extermediate  strain  rates 
and  the  entire  stress-strain  curve  is  obtained  using  strain  gages  and  an  ex- 
tensometer.  The  change  in  diameter  is  used  to  calculate  Poisson's  ratio.  The 
high  strain  rate  tests  are  conducted  on  the  split-Hopkinson  bar  using  the 
standard  tension  test  specimen  for  that  apparatus.  All  of  the  metal  specimens 
have  been  designed  and  fabricated. 

The  advanced  composite  selected  for  the  study  is  boron/aluminum  used  in 
the  APSI  Stage  1  fan  blade  with  a  fiber  layup  of  (0/22/0/-22) .  The  4-mil 
boron  fibers  are  in  a  matrix  of  2024  aluminum.  Each  ply  is  4.7  mils  thick 
and  does  not,  in  general,  run  the  length  of  the  blade.  The  center  ply  of  the 
blade  is  a  plate  of  aluminum.  On  either  side  of  it  is  a  ply  of  stainless 
steel  wire  mesh.  The  next  plies  on  both  sides  are  the  B/Al  layup.  The  entire 
outer  blade  surface  is  covered  with  a  ply  of  stainless  steel  wire  mesh.  The 
leading  edge  has  a  thin  coat  of  pure  nickel. 

The  material  properties  are  needed  in  three  principal  directions  because 
the  composite  is  anisotropic.  The  data  needed  are  the  tensile  stress-strain 
curves  for  each  direction,  three  Poisson's  ratios,  and  three  shear  stress- 
strain  curves.  However,  the  collection  of  this  information  is  simplified  con¬ 
siderably  by  the  use  of  laminate  theory.  This  theory  employs  the  results  of 
material  property  tests  on  unidirectional  specimens  to  calculate  the  bulk 
properties  of  the  angle-ply  layup. 
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The  compos  ice  test  specimens  have  been  designed  and  fabricated  and  Che 
experimental  methods  that  will  be  employed  to  obtain  the  composite  material 
properties  have  been  established. 


6.1.2  Real  Blade  Strain  Rates 


Testing  of  the  full-scale  blades  to  obtain  the  maximum  strain  rates  has 
been  completed.  The  highest  strain  rates  developed  were  less  chan  400  in./ 
in. /sec  in  any  of  the  types  of  blades  tested  (J79,  F101,  and  APSI).  For 
the  J79  impacts,  the  highest  strain  rate  for  the  30%  span  impact  was 
344  in. /in. /sec  for  a  1.5-pound  bird  impact.  This  value  was  for  a  gage 
located  on  the  midchord  at  the  root  in  the  spanwise  direction.  For  the  70% 
span  impact  of  the  J79  blade,  the  highest  strain  rate  was  362  in. /in./ 
sec  which  was  for  a  1.5-pound  bird  impact.  This  value  was  for  a  gage  directly 
behind  the  impact  site  which  was  oriented  in  the  chordwise  direction. 

The  highest  strain  rate  for  the  30%  span  impacts  of  the  tip  restrained 
F101  blade  was  341  in. /in. /sec.  For  the  70%  span  impacts  of  this  blade, 
the  highest  strain  rate  was  367  in. /in. /sec.  Both  strain  rates  were  for 
1.5-pound  bird  Impacts  and  for  a  gage  location  at  the  tip  along  the  leading 
edge  in  the  spanwise  direction. 

For  the  APSI  blade,  the  highest  strain  rate  developed  was  366  in. /in./ 
sec  for  a  70%  span  impact  of  a  3-ounce  bird.  This  occurred  at  a  gage  location 
at  the  midchord  of  the  root  in  the  spanwise  direction.  The  highest  strain  rate 
for  the  30%  span  impact  was  362  in. /in. /sec  for  a  3-ounce  bird.  This  gage 
location  was  behind  the  impact  site  with  the  orientation  of  the  gage  being 
in  the  direction  of  the  chord. 

Figure  9  shows  the  impact  test  setup  used  for  real  blades.  Shown  in  this 
figure  is  the  cantilever  mounted  J79  blade.  The  slicer  and  the  target  blade 
are  oriented  relative  to  the  path  of  the  projectile  to  provide  the  slice 
weight  as  defined  in  Tables  1  and  2.  For  example,  to  simulate  the  impact  on 
the  J79  blade  for  a  rotor  speed  of  7460  rpm  (maximum  engine  speed)  and  a  bird 
velocity  of  200  ft/sec,  bird-blade  incidence  angles  and  impactor  velocities  of 
51.1*  and  560  ft/sec  and  36.4*  and  795  ft/sec  are  required  for  30%  and  70% 
span  impact  locations,  respectively.  For  a  3-ounce  bird,  the  30%  and  70%  span 
impact  parameters  result  in  slice  weights  of  0.126  lb  and  0.129  lb,  respec- 
t ively . 

Figure  10  shows  the  strain  gage  locations  used  to  determine  the  strain 
rates  for  the  J79  blade.  Figures  11  and  12  show  examples  of  the  J79  blade 
response  resulting  from  simulated  bird  impact.  Shown  in  these  figures  are 
the  measured  strain-time  and  computed  strain  rate-time  histories  for  Gages  L 
and  3  for  a  3-ounce  bird  impact  on  the  J79  blade.  Figure  13  shows  the  strain 
in  both  the  time  and  frequency  domains  resulting  from  a  hammer  blow. 
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Figure  9.  Bench  Impact  Test  Setup  for  Real  Blades 
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Strain  and  Strain  Rate  Time  Histories  for  Gage  No.  1 
(at  Root  of  Airfoil  on  Pressure  Side) . 
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Figure  13.  Strain  in  the  Time  and  Frequency  Domain  Resulting  from  Hammer 
Blow  for  Gage  No.  k  J79  Blade,  Test  Date  10/16/78. 


The  bird  impact  response  data  shown  in  Figures  11  and  12  correspond  to 
plastic  response  behavior  in  the  blade.  This  blade  experienced  yielding  at 
both  the  root  and  the  impact  site  as  a  result  of  the  impact.  In  this  case, 
the  actual  impact  velocity  was  984  ft/sec  and  not  795  ft/sec  and  the  actual 
slice  weight  was  0.0619  lb  (28  gms)  and  not  0.129  lb  (58.5  gms)  as  specified 
in  Table  1.  In  addition,  the  substitute  bird  impacted  the  blade  somewhat 
above  the  702  span  location.  The  strain  rates  were  computed  from  the  mea- 
sured  strain-time  history  by  employing  Fourier  transform  techniques.  This 
was  accomplished  by  using  the  duality  between  the  derivative  of  a  time  func¬ 
tion  and  the  complex  product  formed  from  the  Fourier  transform  of  the  time 
function  and  the  factor  (ju)n;  j  is  the  imaginary  direction  in  the  complex 
plane,  u  rad/sec  is  the  frequency  and  n  is  the  order  of  the  derivative.  This 
approach  Is  expressed  mathematically  In  the  development  shown  below. 

The  Fourier  transform  of  the  original  time  function  (strain-time  history) 
X(t)  is  equal  to 


00 

X(u)  -  /  X(c)£“jut  dt 

<*/ 

—  CO 


The  Fourier  transform  Xo(u)  of  the  first  derivative  (strain  rate)  X(t)  can 
be  obtained  from  the  fourier  transform  of  the  original  time  function  X(t)  as 
follows: 


Xq(ui)  ■ 


dt  *  jutX(w) 


Finally,  the  strain  rate,  X(t),  is  obtained  via  the  inverse  Fourier  transform. 


OO 


•  00 


The  spectrum  shown  in  Figure  13  corresponds  to  the  elastic  response  be¬ 
havior  obtained  with  the  hammer  blow.  Two  dominant  response  peaks  are  in  evi¬ 
dence  at  120  Hz  and  350  Hz.  These  frequencies  correspond  to  the  bench  test 
frequencies  for  the  first  and  second  flex  modes  of  the  cantilever-mounted  J79 
blade. 
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6.1.3  Material  Property  Data 


The  material  characterization  tests  on  the  two  metals,  410  stainless 
steel  (410  substituted  for  403)  used  for  the  J79  Stage  1  compressor  blade  and 
titanium  8AI-IM0-IV  used  in  the  F101  Stage  1  fan  blade,  are  currently  under¬ 
way.  The  tests  on  the  advanced  composite  (B/Al)  used  in  the  APSI  Stage  1  fan 
blade  have  not  yet  started.  As  mentioned  previously,  the  advanced  composite 
used  in  the  APSI  blade  is  actually  a  hybrid  composite.  This  blade  has  the 
following  layup:  SSWM/0#/22*/0*/-22*/SSWM/Al/SSWM/-22*/0*/22*/0*/SSWM.  £e 
will,  therefore,  be  necessary  to  characterize  the  stainless  steel  wire  mesh 
also. 


The  strain  rates  for  the  metals  tests  are  shown  in  Table  3.  The  upper 
limit  is  that  value  determined  by  the  real  blade  impact  testing. 


Table  3.  Characterization  of  Metallic  Materials. 


Material 

410  Stainless  Steel 
8AI-IM0-IV  Titanium 


Strain  Rates  (sec  jp 
0.001,  0.1,  1,  10,  100,  Highest 
0.1,  10,  Highest 


The  metal  tests  at  the  lower  strain  rates,  0.001  e/sec,  0.1  e/sec,  and 
1.0  e/sec,  have  been  conducted.  In  addition,  the  high  strain  rate  material 
properties  tests  using  the  split  Hopkinson  bar  apparatus  have  also  been  con¬ 
ducted.  The  strain  rates  for  these  high  strain  rate  tests  varied  from  330  to 
570  e/sec  for  the  8-1-1  titanium,,  and  340  to  760  e/sec  for  the  stainless 
steel  (the  larger  range  in  strain  rates  on  the  stainless  steel  was  achieved 
by  setting  the  striker  bar  to  attain  different  velocities).  The  maximum  rates 
are  higher  than  the  impact  tests  on  the  blades  indicated  were  necessary.  The 
data  have  been  digitized  and  graphed.  The  metal  tests  at  10  e/sec  and 
100  e/sec  have  not  yet  been  accomplished.  The  low  and  intermediate  strain 
rate  tests  (0.001-10  e/sec)  have  been  conducted  per  standard  ASTM  E8-69?  and 
the  specimen  definition  is  shown  in  Figure  14.  The  standard  tensile  specimen 
used  for  the  high  strain  rate  tests  on  the  split  Hopkinson  bar  is  shown  in  Fig¬ 
ure  15. 


The  test  matrix  that  will  be  used  for  the  characterization  of  the  boron- 
aluminum  and  the  stainless  steel  wire  mesh  is  shown  in  Table  4.  For  the  ten¬ 
sile  tests  in  the  direction  of  the  symmetry  axes  (0*  and  90*),  a  modified 
ITTRI  specimen  configuration  will  be  used.  This  specimen  is  shown  in  Fig¬ 
ure  16.  For  the  shear  tests  in  the  interlaminar  plane,  the  10°  off-axis  shear 
test  developed  by  Chamis  and  Sinclair^  will  be  used  for  the  B/Al.  A  drawing 
of  the  specimen  appears  in  Figure  17.  Torsion  rod  specimens  will  be  used  to 
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250  in.  x  0.005 


Figure  15.  Split  Hopkinson  Bar  Specimen. 


Figure  16.  Modified  ITTRI  Specimen. 


determine  the  transverse  shear  properties  for  both  the  B/Al  and  the  stainless 
steel  wire  mesh.  The  in-plane  shear  properties  for  the  stainless  steel  wire 
mesh  will  be  determined  with  the  ±45*  shear  panel  specimen. 


Table  4.  Characterization  Tests  for  B/Al  and 
Stainless  Steel  Wire  Mesh. 


Material 

Test  Description 

Strain  Rates 

B/Al 

0*  Tension 

1,10 

90*  Tension 

1,10 

10*  Off-axis  Shear 

10 

Rod  Torsion  Shear 

10 

SSWM 

Tension  (in  plane) 
Torsion  (through  the 

1,10 

thickness) 

1,10 

±45  Shear 

10 

Rod  Torsion  Shear 

10 

The  through-the-thickness  tensile  specimen  shown  in  Figure  18  will  be 
used  to  determine  the  normal  to  the  plane  tensile  properties  for  the  stain¬ 
less  steel  wire  mesh.  The  tensile  specimens  (include  10*  off-axis  and  ±45* 
panels)  will  be  instrumented  with  high  resistance  strain  gages.  For  the 
through-the-thickness  and  torsion  rod  specimens,  it  is  planned  to  use  an  ex- 
tensometer  to  measure  displacement. 

It  is  not  anticipated  that  the  composite  will  be  significantly  strain 
rate  sensitive.  Strain  rate  effects  in  B/Al  were  investigated  by  Krinke, 
Barber,  and  Nicholas^  using  a  three-point  bend  (Charpy  impact)  technique. 
They  found  no  significant  strain  rate  effects  over  6  orders  of  magnitude;  the 
highest  strain  rate  employed  in  the  study  was  100  sec-*. 


6.2  SUBTASK  B  -  LOCAL  LEADING  EDGE  DAMAGE 
6.2.1  Introduction 

This  subtask  is  concerned  with  the  generation  of  empirical  design  data 
for  local  leading  edge  damage  to  first  stage  fan  and  compressor  blades. 
Testing  is  being  conducted  to  determine  and  quantify  the  damage  caused  by 
leading  edge  impacts  for  a  range  of  pertinent  impact  conditions.  The  results 
of  this  study  will  include  data  expressed  in  terms  of  the  residual  properties 
of  the  materials  which  are  directly  relevant  to  fan  and  compressor  designers. 
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6.2.2  Overall  Approach 


Three  types  of  blade  materials,  geometries,  and  sizes  will  be  investi¬ 
gated  using  substitute  bird  materials,  ice,  granite  pebbles  and  steel  spheres 
as  the  impactors.  This  range  of  impactors  is  representative  of  the  types  of 
objects  commonly  ingested  into  engines.  The  three  blade  types  are  the  J79 
blade  using  type  403  stainless  steel,  the  F101  blade  using  8-1-1  titanium, 
and  the  APSI  metal  matrix  boron/ aluminum  blade.  Fabricated  metal  and  compos¬ 
ite  specimens  and  actual  APSI  blades  will  be  tested.  The  geometries  of  the 
leading  edges  of  the  fabricated  specimens  will  be  similar  to  the  geometries 
of  the  blade  types  at  the  50%  span  location.  Three  specific  investi¬ 
gations  are  planned  in  this  subtask.  The  first  establishes  a  base  of  damage 
data.  The  second  addresses  scaling  and  subscale  tests,  while  the  third  in¬ 
vestigates  the  effects  of  preload  and  possible  damage  enhancement.  In  con¬ 
junction  with  the  above  investigations,  fatigue  and  tensile  tests  will  be  per¬ 
formed  on  damaged  and  undamaged  specimens  to  establish  the  effect  of  impact 
on  these  properties. 

Damage  base  data  will  be  generated  using  nominal  thickness  fabricated 
specimens  (the  leading  edge  thickness  will  be  identical  to  the  corresponding 
blade  type  at  the  50%  span  level).  Half -thickness  fabricated  specimens 
will  then  be  used  in  a  scaling  effects  investigation.  Preload  investigations 
will  be  performed  on  the  fabricated  metal  specimens  to  simulate  the  effects 
of  centrifugal  loading  on  local  leading  edge  damage.  The  specimens  will  be 
preloaded  at  stresses  typical  of  the  impact  locations  (30  and  70%  span  levels). 


6.2.3  Metal  Specimen  Tests 

The  impactors  for  the  metal  specimen  tests  include  1.5-pound  birds, 
1-inch  diameter  ice  spheres,  1-inch  diameter  by  4-inch  long  ice  cylinders, 
0.062  and  0.125-inch  diameter  steel  spheres  and  1/4-inch  diameter  pebbles. 

The  impact  velocities  that  will  be  used  for  the  metal  specimens  corre¬ 
spond  to  those  which  would  be  typical  of  an  impact  at  the  70%  span  and  at  the 
30%  span  at  full  power  settings  of  the  engine  during  takeoff  on  each  of  the 
metal  blade  types  (refer  to  Table  1).  Impacts  at  70%  span  are  representative 
of  the  highest  velocity  impacts  experienced  by  a  blade.  Impacts  at  30%  span 
are  typical  of  those  in  the  highest  stress  regions  of  the  blade  where  it  is 
most  vulnerable  to  the  effects  of  impact  degradation. 


The  hard  body  impact  tests  on  the  nominal  and  half-scale  thickness  strain¬ 
less  steel  fabricated  specimens  have  been  completed.  The  soft-body  projectile 
impact  testing  using  3-ounce  microballoon  gelatin  birds  was  initiated  for  the 
stainless  steel  specimens. 
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6.2.4  Composite  Specimen  Tests 

6. 2. 4.1  Laboratory  Size  Composite  Specimen  Tests 

Fabricated  B/Al  composite  specimens  that  simulate  the  leading  edge  of 
the  APSI  blade  at  50%  span  will  be  tested.  One  impact  velocity  that  corre¬ 
sponds  to  that  typical  of  an  impact  at  the  70%  span  location  of  the  APSI 
B/Al  blade  at  the  takeoff  power  setting  will  be  used  (refer  to  Table  1) . 

In  the  damage  data  base  study,  data  will  be  obtained  for  two  sizes  (0.062 
and  0.125-inch  diameter)  of  steel  sphere  impactors.  In  the  scale  effects 
investigation,  data  will  be  obtained  only  with  the  0.062-inch  steel  sphere 
impactor.  No  testing  is  planned  in  regard  to  preload  effects,  since 
insignificant  effect  of  preload  was  experienced  for  center  panel  impacts  onto 
composite  specimens  in  a  previous  investigation.  Seven  undamaged  specimen 
tests  will  be  used  to  provide  a  data  base  for  the  residual  property  measure¬ 
ments.  A  total  of  22  fabricated  B/Al  composite  specimens  will  be  used. 


6. 2. 4. 2  Full-Scale  Composite  Blade  Tests 

In  addition  to  testing  laboratory  size  specimens,  full-scale  composite 
blades  will  also  be  impact  tested.  The  local  leading  edge  damage  problem  will 
be  investigated  by  conducting  hard  impactor  impacts  on  the  leading  edge  of 
three  B/Al  APSI  blades.  The  hard  impactors  utilized  in  the  testing  will  be 
0.25-inch  diameter  pebbles,  0.062-inch  diameter  steel  spheres,  and  0.125-inch 
steel  spheres.  Mulitiple  impacts  along  the  leading  edge  of  the  real  blades 
will  allow  determination  of  the  ballistic  limit  and  provide  qualitative 
damage  data  for  the  three  sizes  of  hard  impactors. 

Full-scale  composite  blade  testing  will  also  be  conducted  to  investigate 
structural  damage  resulting  from  soft  body  projectile  impacts.  Four  B/Al  APSI 
blades  will  be  impacted  using  either  3-ounce  microballoon  gelatin  birds  or 
1-inch  diameter  ice  spheres.  The  blades  will  be  impacted  at  the  30  and  70% 
span  locations  at  velocities  and  impact  angles  corresponding  to  the  span 
location. 


6.3  SUBTASK  C  -  MATERIAL  SCREENING  TESTS 
6.3.1  Introduction 


This  subtask  addresses  the  development  of  material  screening  tests  neces¬ 
sary  to  evaluate  and  effectively  rank  candidate  materials  for  fan  and  com¬ 
pressor  blades  that  possess  superior  F0D  resistance.  The  material  parameters 
that  will  be  investigated  in  this  study  include:  (1)  the  perforation  resis¬ 
tance  of  the  material;  (2)  the  extent  the  material  may  plastically  deform; 

(3)  the  extent  to  which  the  material  is  vulnerable  to  catastrophic  structural 
failure;  and  (4)  the  extent  to  which  the  material  is  vulnerable  to  degradation 
of  fatigue  properties. 
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In  this  investigation,  it  is  planned  to  define  these  material  parameters 
in  terms  of  response  to  impact  in  small-scale  tests.  The  perforation  resis¬ 
tance  material  parameter  will  be  quantified  as  a  ballistic  limit  velocity,  Vl» 
for  a  given  impact  configuration.  The  quantity  is  a  direct  measure  of  the 
amount  of  impact  energy  which  the  material  can  absorb  without  catastrophic 
local  damage.  The  plastic  deformation  parameter  will  be  quantified  by  con¬ 
ducting  impact  tests  which  maximize  the  likelihood  of  plastic  deformation. 

The  catastrophic  structural  failure  material  parameter  will  be  quantified  in 
a  Charpy  impact  test  designed  to  produce  gross  structural  failure.  The  fa¬ 
tigue  properties  parameter  will  be  quantified  in  terms  of  reduction  of  ulti¬ 
mate  fatigue  strength.  The  parameters  defined  by  these  tests  may  be  con¬ 
sidered  individually  or  combined  with  suitable  weighting  functions  to  compute 
an  overall  figure  of  merit  for  a  particular  material. 


6.3.2  Materials  to  be  Investigated 

The  usefulness  of  material  ranking  by  the  screening  tests  will  be  eval¬ 
uated  by  applying  the  tests  to  six  materials.  Three  of  the  materials  will  be 
studied  extensively  in  the  overall  program.  These  include  403  stainless  steel 
used  in  J79  Stage  1  compressor  blades,  8-1-1  titanium  used  in  F101  Stage  l  fan 
blades,  and  boron/ aluminum  used  in  APSI  Stage  1  fan  blade.  Additional  mate¬ 
rials  to  be  tested  include  a  heat  treated  403  stainless  steel,  6A1-4V  titanium, 
and  a  graphite  epoxy  composite.  All  of  the  specimens  (metal  and  composite) 
have  been  fabricated. 


6.3.3  Screening  Tests 

6. 3. 3.1  Ballistic  Limit  Testing 

This  test  will  evaluate  the  ballistic  limit  velocity  (Vl>  for  normal 
impacts  of  various  typical  FOD  projectiles.  Three  impactors  will  be  used  in 
the  testing.  A  1/2-inch  diameter  microballoon  gelatin  sphere  with  10  to 
15%  porosity  will  be  utilized  on  all  six  materials  for  Bird  simulation. 

The  remaining  two  impactors  are  1/2-inch  ice  balls  and  i/4-inch  diameter 
stones  or  "standard"  pebble.  Both  ice  and  stone  impacts  will  be  used  on  two 
of  the  metals  which  will  be  studied  extensively  in  the  overall  program. 

These  will  be  the  403  stainless  steel  and  the  8-1-1  titanium.  Pebble  impacts 
will  also  be  conducted  on  the  boron/aluminum  composite  specimens;  however, 
ice  impacts  will  not  be  conducted  since  ice  is  believed  to  behave  similar  to 
the  substitute  bird  material  (microballoon  gelatin).  A  total  of  25  composite 
specimens  will  be  used  for  ballistic  limit  testing.  The  launch  velocity  of 
the  impactors  will  te  varied  until  the  ballistic  limit  (Vl)  is  determined 
within  about  5%. 

A  number  of  ballistic  limit  tests  (using  the  granite  pebble)  has  been 
conducted  on  the  metal  and  composite  specimens.  The  ballistic  limit  of  the 
403  stainless  steel  was  determined  to  be  about  3375  feet  per  second  for  the 
1/4-inch  diameter  granit  pebbles.  The  ballistic  limit  for  pebble  impicts  on 
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Che  boron/aluminum  specimens  was  determined  Co  be  about  620  feec  per  second. 
Thus,  Che  ballistic  limic  of  scainless  steel  is  five  times  greater  than  for 
boron/aluminum  for  peb  le  impacts. 


6. 3. 3. 2  Maximum  Deformation  Testing 

Specimens  that  are  identical  to  those  used  in  the  ballistic  limit  test 
will  be  utilized  in  the  maximim  deformation  testing.  The  specimens  will  be 
impacted  by  1/2-inch  diameter  microballoon  gelatin  birds  (substitute  bird 
material)  to  maximize  plastic  deformation  of  the  target.  Again,  normal  im¬ 
pacts  in  the  specimen  center  will  be  utilized.  The  velocities  will  also  be 
varied  to  collect  data  such  that  a  plot  of  the  deformation  damage  versus 
velocity  can  be  made.  It  is  believed  that  the  maximum  deformation  damage  will 
occur  at  a  velocity  just  below  the  ballistic  limit  velocity  (V^),  probably 
0.9  VL. 


The  results  of  the  local  leading  edge  damage  effort  should  provide  a 
basis  for  selecting  the  most  diagnostic  damage  parameter  to  use  in  measuring 
the  damage.  TWo  likely  candidates  are  the  plastically  deformed  area  and  max¬ 
imum  deflection.  A  total  of  10  composite  specimens  will  be  used  for  maximum 
deformation  testing. 

6. 3. 3. 3  Gross  Structural  Damage  Testing 

This  test  is  designed  to  assess  vulnerability  to  catastrophic  structural 
damage  of  blades  by  using  test  specimens  from  the  six  candidate  materials. 

Such  damage  is  probably  due  to  bending  stresses  at  the  root  or  some  other  mode 
on  an  impacted  blade.  The  results  of  failure  mode  tests  from  previous  tasks 
of  the  overall  program  will  clarify  this  point.  It  is  anticipated  that  an  in¬ 
strumented  Charpy  impact  test  will  be  sufficient  to  characterize  material  re¬ 
sponse  to  this  type  of  loading.  A  total  of  eight  composite  specimens  will  be  used. 


6. 3. 3. 4  Fatigue  Testing 

The  fatigue  strength  test  is  an  attempt  to  classify  damage  to  specimens 
in  terms  of  an  equivalent  stress  concentration  factor.  Holes  and  cracks  will 
be  machined  in  specimens  using  the  electrical  discharge  machining  technique 
(EDM).  The  holes  (0.177-in.  diameter)  and  cracks  (0.177-in.  long)  will  be  machined 
in  the  center  of  the  specimen  at  its  midspan.  The  damaged  specimens  will  then 
be  tested  in  a  tensile  fatigue  machine.  Each  group  of  specimens  with  identi¬ 
cal  damage  will  be  tested  in  the  fatigue  machine  at  load  levels  to  obtain 
failure  in  the  range  from  10^  to  10^  cycles.  The  number  of  cycles  to  failure 
(complete  separation)  will  be  recorded  for  each  specimen.  The  data  will  be 
plotted  in  the  form  of  net  section  average  stress  against  number  of  cycles  to 
failure  for  each  group  of  specimens.  One  or  two  specimens  of  each  group  will 
be  used  to  obtain  the  residual  tensile  strength  for  baseline  purposes.  A 
total  of  24  composite  specimens  will  be  used. 
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7.0  TASK  V  -  PARAMETRIC  ANALYSIS 


7 . 1  INTRODUCTION 


Parametric  relationships  describing  the  changes  in  dynamic  structural 
response  of  impacted  simple  plate  elements  with  the  progressive  introduction 
of  blade  airfoil  geometric  features  will  be  derived  in  this  task.  Addition*- 
ally,  the  adequacy  of  the  structural  response  models  formulated  in  Tasks  II, 
III,  and  IV  will  be  investigated  and  modifications  made  as  necessary. 

This  task  is  broadly  divided  into  three  phases.  The  first  phase  is  anal¬ 
ysis  and  correlation  of  Task  VI  tests  conducted  with  simple  cantilevered  flat 
plates  with  different  aspect  ratios  and  thicknesses.  The  second  phase  is  the 
analysis  and  correlation  of  Task  VI  tests  with  geometry  and  impact  conditions 
which  are  more  representative  of  actual  airfoil  FOD  dynamics  but  which  are 
still  simple  enough  that  the  sources  of  observed  effects  can  be  isolated. 

The  third  and  final  phase  of  this  task  is  the  analysis  and  correlation  of 
Task  VI  and  Task  IX  tests  of  the  three  selected  first-stage  airfoils  in  sta¬ 
tionary  and  rotating  environments. 


7.2  BLADES  SELECTED  FOR  ANALYSIS 


The  three  first-stage  blades  selected  are  representative  of  first-stage 
airfoils  from  Air  Force  inventory,  production  development,  and  advanced  de¬ 
velopment  engines.  These  blades  include  subsonic,  transonic,  and  supersonic 
airfoil  shapes  and  stainless  steel,  titanium,  and  an  advanced  composite  mate¬ 
rial.  The  blades  selected  are  listed  below: 


Blade 

Representing 

Airfoil 

Material 

Shrouded 

J79 

Air  Force  Inventory 

Subsonic 

Stainless  Steel 

No 

F101 

Production  Development 

Transonic 

Titanium 

Yes 

APSI 

Advanced  Development 
Supersonic  Engines 

Supersonic 

B/Al  Composite 

No 

7.3  PARAMETRIC  MATRIX 


A  parametric  matrix  has  been  formulated  to  define  the  impact  conditions 
and  response  model  geometries  that  will  be  analyzed  in  Tasks  V  and  X  and  pro¬ 
vide  guidance  in  the  selection  of  the  structural  element  and  full-scale  blade 
tests  that  will  be  conducted  in  Tasks  VI,  IX  and  XI.  This  matrix  is  shown  in 
Table  5. 

Table  5  describes  the  structural  elements,  element  fixity  and  material, 
the  analysis  level,  the  loading  and  impactor,  and  the  impact  location  and 
incidence.  The  analysis  level  in  Table  5  refers  to  the  transient  structural 
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Table  5.  Parametric  Matrix  Defining  Structural  Elements  and  Impact 
Conditions  to  be  Analyzed. 
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Table  5.  Parametric  Matrix  Defining  Structural  Elements  and  Impact 
Conditions  to  be  Analyzed  (Continued). 
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Table  5.  Parametric  Matrix  Defining  Structural  Elements  and  Impact 
Conditions  to  be  Analyzed  (Concluded) . 


Croup 

Structural  Clement 
and  Comments 

Analysis 

Level 

Loading/ 

Inpactor 

Inpactor  Location 

Inpact  Shroud 

Incidence  Restraint 

Specimen 

Material 

Task 

JO 
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Nat.  fren.-mode  shapes 

2 

HA 

HA 

KA 

NO 
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HA 
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n 

63 
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\i 
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NA 

response  computer  code;  Level  1  Is  the  finite  element  model  based  on  the 
NOSAPM  program,  and  Level  2  Is  the  component  element  model  based  on  the 
COMET  program.  Analyses  of  the  same  structural  model  with  the  Level  1  and 
Level  2  analyses  are  Indicated  In  Table  5  for  a  number  of  cases.  If  the 
agreement  between  the  Level  1  and  Level  2  analyses  proves  to  be  excep¬ 
tionally  good,  then  the  number  of  Level  1  analyses  may  be  reduced. 

Table  5  includes  a  description  of  the  structural  models  and  impact  con¬ 
ditions  that  will  be  used  in  Task  X.  In  Task  X,  Full  Stage  Response  Predic¬ 
tions,  the  response  of  a  full  stage  of  blades  will  be  predicted.  In  this 
case,  it  is  assumed  that  the  only  mechanical  coupling  is  through  shroud-to- 
shroud  interaction  and  the  models  will  include  a  single  blade  with  the  shrouds 
represented  by  springs  and  a  sector  of  shrouded  blades.  The  blades  used  in 
these  models  will  be  the  F101  tip  shrouded  titanium  first  stage  fan  blades 
rotating  at  the  take  off  angular  velocity  of  7550  rpm. 

Details  of  the  shape,  size  and  configuration  of  the  structural  elements 
are  provided  in  Table  6. 


7.4  MODELING  OF  THE  STRUCTURAL  SPECIMENS 


A  study  has  been  initiated  to  consider  several  of  the  options  open  to 
the  user  of  the  first  level  response  analysis  code  (NOSAPM)  which  have  a  bear¬ 
ing  on  the  accuracy  of  the  program's  results.  One  of  the  options  available  to 
the  user  is  the  number  of  nodes  he  chooses  tu  describe  the  geometry  and  the 
displacement  filed  of  the  three  dimensional  fini  e  elements.  A  maximum  of  21 
nodes  per  element  are  available.  Professor  Bathe,  the  author  of  ADINA  and 
consultant  on  this  program,  has  suggested  that  16  nudes  be  used,  eliminating 
the  mid-side  nodes  through  the  thickness.  An  initial  comparison  of  these  two 
element  descriptions  in  the  case  of  a  static  bending  situation  is  made  in  Fig¬ 
ures  19  and  20.  The  problem  used  for  comparison  is  a  cantilever  plate  loaded  at 
its  free  end  with  a  total  force  of  60  pounds  distributed  uniformly  across  its 
total  width.  The  dimensions  of  the  plate  are  7  inches  by  3  inches  by  0.1875 
inch,  as  shown  in  Figure  19.  Four  20-noded  elements  were  used  in  the  finite 
element  analysis  of  this  problem,  and  the  resulting  bending  stresses  incurred 
at  the  top  surface  of  the  plate  both  along  the  centerline  and  at  the  y  *  1.5 
edge  are  plotted  in  Figure  19  as  a  function  of  z  and  compared  with  the  beam  anal¬ 
ysis  solution.  The  solutions  are  in  good  agreement  except  in  the  vicinity  of 
the  clamped  end  where  the  finite  element  analysis  predicts  an  end  effect.  A 
second  finite  element  analysis  was  then  carried  out  using  16-noded  instead  of 
20-noded  elements,  and  the  stress  results  of  this  analysis  are  compared 
against  those  obtained  from  the  20-noded  element  analysis  in  Figure  20.  As 
can  be  seen  from  the  figure,  the  agreement  is  quite  close  except  in  the  vicin¬ 
ity  of  the  clamped  edge.  Additional  comparisons  of  these  element  descriptions 
will  be  subsequently  carried  out.  Since  elimination  of  four  mid-side  nodes 
through  the  thickness  of  each  element  will  lead  to  substantial  cost  savings, 
the  use  of  16-noded  elements  should  be  closely  considered. 


51 


Table  6.  Shape,  Size  and  Configuration  Details  of  Structural  Elements. 
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Figure  20.  Comparison  of  16-Noded  Finite  Element  Solution  with  20-Noded 
Element  Solution. 
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In  addition  to  this  comparison,  another  study  has  been  initiated  to 
determine  the  number  of  elements  that  will  be  necessary  to  adequately  model 
the  plates  and  blades  that  will  be  tested  in  this  program.  The  normal  impact 
of  a  0.0984-inch  thick  titanium  plate  (half  the  thickness  representative  of 
the  F101  blade)  is  being  considered.  Various  impact  velocities  of  a 
1-1/2-inch  diameter,  3-ounce  bird  are  being  considered  and  subsequent  blade 
responses  studied. 
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8.0  TASK  VI  -  STRUCTURAL  ELEMENT  TESTS 


8 . I  INTRODUCTION 


The  work  of  Task  VI,  "Structural  Element  Tests",  is  concerned  with  per¬ 
forming  nonrotating  bench  impact  tests  on  test  specimens  ranging  from  simple 
elements  to  real  blades.  The  simple  elements,  such  as  beams  and  plates,  will 
be  tested  with  progressive  introduction  of  airfoil  geometric  parameters  to 
validate  experimentally  the  analytical  predictions  of  Task  V  and  VIII.  The 
purpose  of  Task  V  is  to  derive  parametric  relationships  describing  the  changes 
in  dynamic  structural  response  as  a  function  of  material  and  geometric  fea¬ 
tures.  The  purpose  of  Task  VIII  is  to  derive  criteria  for  predicting  blade 
foreign  object  impact  damage  tolerance.  These  criteria  will  be  formulated  to 
make  full  use  of  the  transient  dynamic  analysis  and  experimental  results  ob¬ 
tained  from  other  tasks  of  the  program. 


8.2  COMPOSITE  SPECIMEN  TESTS 


8.2.1  Laboratory  Size  Composite  Specimen  Tests 

The  total  number  of  fabricated  B/Al  specimens  that  will  be  tested  is 
equal  to  48.  These  specimens  include  flat  panel,  airfoil-like  cross  section, 
cambered,  and  cambered  and  twisted  shapes.  The  test  conditions  are  listed  in 
Table  7  and  the  B/Al  composite  specimen  configurations  are  shown  in  Figure  21. 


Table  7.  Task  VI  Fabricated  Composite  Specimen  Test  Conditions. 


Velocit  tes : 

Impact  Angle: 

Impactor  2  to  3-Ounce  Bird: 


No  Damage 
Threshold  Damage 
Severe  Damage 

702  Span  Location  -  Oblique 
302  Span  Location  -  Oblique 

All  Configurations  Except 


Configuration  4 

Test  Conditions  for  2  to  3-Ounce  Bird  Impacts  -  36 
Impactor  2-Inch  Ice  Ball:  Configurations  6  and  7 

Test  Conditions  for  Ice  Ball  Impacts  -  12 

Total  Test  Conditions  -  48 
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Figure  21.  Task  VI  B/Al  Composite  Test  Specimen  Coni  iuurat  i>m.- . 
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8.2.2  Full-Scale  Composite  Blade  Tests 

In  addition  to  testing  of  laboratory  size  composite  specimens,  a  number 
of  tests  will  be  conducted  on  the  full  scale  B/Al  APSI  blades.  The  series  of 
tests  to  be  performed  on  real  composite  blades  is  outlined  in  Table  8. 


Table  8.  Task  VI  -  Real  Composite  Blade  Test  Conditions. 


Impact  Parameters:  Mass  -  2  to  3-Ounce  Bird 

-  2-Inch  Ice  Ball 

Velocity  -  No  Damage 

-  Threshold  Damage 

-  Severe  Damage 

Location/ Angle  -  70%  Span/ Edge-Oblique  .  . 

-  30%  Span/ Edge-Oblique 

Test  Conditions  for  Composite  Blades  (12) 


(2) 

(3) 


3.3  METAL  SPECIMEN  TESTS 


Fabricated  titanium  8A1-1MO-1V  and  410  stainless  steel  specimens  and 
the  corresponding  blade  types  (F101  Stage  1  fan  blade  and  J79  Stage  1  compres¬ 
sor  blade)  will  be  tested. 

A  baseline  series  of  tests  will  be  conducted  on  the  titanium  material 
and  a  supplementary  series  of  tests  will  be  conducted  on  the  stainless  steel. 
Titanium  is  chosen  as  the  baseline  material  as  it  is  the  most  common  current 
blade  material. 

The  baseline  series  of  tests  on  the  titanium  material  will  consider  all 
impact  conditions  and  blade  geometrical  effects  introduced  progressively  (not 
independently)  except  camber  and  twist.  It  has  been  established  that  camber 
and  twist  are  not  cold-formable  on  the  titanium  specimens  and  that  these  pa¬ 
rameters  could  be  incorporated  only  through  hoc  forming  or  machining  opera¬ 
tions.  Therefore,  camber  and  twist  will  be  investigated  utilizing  stainless 
steel  specimens  which  will  be  cold  formed. 

Table  9  shows  the  test  conditions  for  fabricated  metal  specimens  and 
Table  10  shows  the  test  conditions  for  Che  real  metal  blades. 

During  testing  of  the  fabricated  structural  specimens,  the  different 
geometry  effects  (i.e.,  aspect  ratio,  thickness/chord  ratio,  etc.)  wilL  be 
introduced  in  a  different  sequence  for  the  titanium  and  steel  specimens  so 
that  any  synergistic  effects  that  might  exist  can  be  accurately  assessec. 
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Table  9.  Task  VI  -  Fabricated  Metal  Specimen  Test  Conditions. 


Baseline  Test  Conditions  for  Titanium 


Impact  Parameters: 


Blade  Parameters: 


Mass  (All)  (3) 

Velocity  (All)  (3) 

Location/Angle  (All)  (2) 

All  (Progressive)  Except  Camber  (8) 

and  Twist.  Also  One  Parameter 
at  Elevated  Temperature.  Also 
Center  and  Edge  Impact  on  Blade 
Type  Aspect  Ratio/Panels. 


Test  Conditions  -  144 


Supplementary  Series  for  Stainless  Steel 


Mass  -  1.5-Pound  Bird 

-  2  to  3-Ounce  Bird  (2) 

Impact  Parameters:  Velocity  -  (All)  (3) 

Location/Angle  -  (All)  (2) 

Blade  Parameters:  Three  Most  Sensitive  from  (5) 

Experiments  on  Titanium 
Plus  Camber  and  Twist 


Test  Conditions  -  60 


Table  10.  Task  VI  -  Real  Metal  Blade  Test  Conditions. 


Impact  Parameters: 

Mass  -  2  to  3-Ounce  Bird 

(3) 

-  1-1 /2-Pound  Bird 

-  2-Inch  Ice  Ball 

Velocity  -  All 

(3) 

Location/ Angle  -  All 

(2) 

Blades: 

F101  (Shrouded) 

J79 

(2)-Metal 

Blades 

Test  Conditions  for 

Metal  Blades 

(36) 

9.0  TASK  VII  -  ERROR  BAND  ANALYSIS 


9.1  INTRODUCTION 

The  objective  of  this  task  is  to  analytically  determine  the  uncertain* 
ties  in  strain  history,  deflection  history,  and  damage-related  parameters  of 
FOD  which  are  introduced  by  such  factors  as  manufacturing  tolerance,  variabil¬ 
ity  in  material  properties,  and  variability  in  impact  conditions.  This  deter¬ 
mination  is  important  in  Task  VIII  for  quantifying  variances  for  the  damage 
model  and  in  Tasks  VI  and  IX  for  comparing  computer  output  with  test  results. 

The  damage-related  parameters  will  be  among  those  used  in  Task  VIII  as 
the  basis  for  the  impact  damage  model,  for  example,  maximum  stress  or  maximum 
strain.  These  parameters  will  be  obtained  as  output  from  the  analysis  models 
defined  in  Task  II.  It  is  expected  that  the  second  level  response  model  will 
be  adequate,  but  the  detailed  modeling  required  for  studying  the  effects  of 
local  deviations  in  geometry  within  manufacturing  tolerance  might  require  the 
use  of  the  first  level  analysis  response  model. 


9.2  PROPERTIES  REQUIRING  AN  ERROR  BAND  ESTIMATION 

A  list  of  properties  relating  to  blade  material,  impactor  characteris¬ 
tics,  and  aerogeometric  factors  for  which  the  FOD  design  criteria  developed 
in  Task  VIII  will  most  need  an  error  band  estimation  has  been  compiled. 

It  is  believed  that  variation  of  these  properties  within  a  reasonable  range 
can  cause  a  wide  variation  in  damage  resulting  from  foreign  object  impact. 
These  error  band  properties  and  the  associated  damage  modes  are  identified 
in  Che  matrix  shown  on  Table  11. 


Table  11.  Matrix  of  Error-Band  Properties  and  Associated  Damage  Modes 


10.0  TASK  VIII  -  FOREIGN  OBJECTIVE  IMPACT  DESIGN  CRITERIA 


10.1  INTRODUCTION 


The  purpose  of  this  cask  is  Co  derive  Che  criCeria  for  prediccing  blade 
foreign  objecc  impacC  damage  Colerance.  The  criCeria  will  be  formulaced  Co 
make  full  use  of  Che  CransienC  dynamic  analysis  and  experimental  resulCs 
obcained  in  Che  previous  casks  and  will  be  based  on  a  comb ina cion  of  cheory 
and  regression  analysis.  Dimensional  analysis  and  an  examinacion  of  Che 
physics  of  impacC  will  decermine  which  impacC  dynamics  and  airfoil  response 
paramecers  will  be  used  Co  formulaCe  Che  damage  model.  Regression  analysis 
will  be  used  co  decermine  Che  consCanCs  in  Che  model  by  correlacion  with  Che 
experimencal  resulCs.  The  experience  gained  at  General  Electric  in  the 
development  of  simplified  criCeria  based  on  purely  elasCic  analyses  will  pro¬ 
vide  an  excellent  basis  for  Che  complete  elastic-plastic  analysis  of  chis 
program. 

The  damage  models  formulaced  will  be  correlated  with  Che  observed  damage 
from  Tasks  VI  and  IX  tesCs,  and  will  each  be  sufficiencly  resolved  Co  predict 
a  chreshold  value  beneach  which  Che  damage  is  noc  observed  and  above  which 
the  prescribed  damage  is  evident.  The  level  of  damage  for  such  definitions 
as  local  plastic  deformation  will  be  of  relatively  increasing  severity  from 
no  damage  Co  large  plascic  deformation  so  ChaC  multiple  qualitative  thresholds 
will  be  required.  A  more  definitive  damage  criterion  is  to  specify  Che  size 
of  the  region  of  local  plastic  deformation  relative  Co  a  normalizing  dimension 
of  Che  airfoil  such  as  Che  chord.  Damage  definicions  such  as  local  fracCure 
are  qualitatively  more  precise  measures  so  thac  a  single  chreshold  will  be 
sufficient . 

10.2  FOREIGN  OBJECT  IMPACT  DESIGN  CRITERIA  GOALS 


An  inicial  list  of  foreign  object  impact  design  criteria  goals  has  been 
compiled  for  consideration.  For  both  metal  and  composite  blades,  these  cri¬ 
teria  include  limits  on  overall  blade  displacement  to  avoid  interference, 
limit  on  local  plastic  deformation  to  avoid  stall/surge,  and  limits  on  the 
total  distortion  strain  energy  per  unit  volume  to  avoid  mass  loss.  Addi¬ 
tional  criteria  for  composite  blades  include  limits  on  the  shear  and  through- 
the-thickness  normal  stresses  to  avoid  delamination,  and  limits  on  the  span- 
wise  and  chordwise  normal  stresses  to  avoid  fiber  failure. 

Relative  to  metal  blades,  we  are  looking  at  an  approach  which  appears 
very  promising  for  predicting  tearing  that  is  based  on  the  relative  values 
of  the  maximum  and  minimum  principal  strains.  This  approach  is  based  on  a 
concept  used  in  forming  sheet  metal  where  a  forming  limit  diagram  is  used 
to  determine  whether  a  metal  can  be  shaped  without  risk  of  failure  or  tear¬ 
ing. 
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